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ABSTRACT 

We carry out an independent search of Kepler photometry for small transiting planets with sizes 
0.5-8.0 times that of Earth and orbital periods between 5 and 50 days, with the goal of measuring 
the fraction of stars harboring such planets. We use a new transit search algorithm, TERRA, optimized 
to detect small planets around photometrically quiet stars. We restrict our stellar sample to include 
the 12,000 stars having the lowest photometric noise in the Kepler survey, thereby maximizing the 
detectability of Earth-size planets. We report 129 planet candidates having radii less than 6 Re 
found in 3 years oi Kepler photometry (quarters 1-12). Forty-seven of these candidates are not in 
iBatalha et al.l (|2012l ), which only analyzed photometry from quarters 1-6. We gather Keck HIRES 
spectra for the majority of these targets leading to precise stellar radii and hence precise planet 
radii. We make a detailed measurement of the completeness of our planet search. We inject synthetic 
dimmings from mock transiting planets into the actual Kepler photometry. We then analyze that 
injected photometry with our TERRA pipeline to assess our detection completeness for planets of 
different sizes and orbital periods. We compute the occurrence of planets as a function of planet 
radius and period, correcting for the detection completeness as well as the geometric probability of 
transit, R^/a. The resulting distribut ion of planet s i zes ex hibits a power law rise in occurrence from 
5.7 Re down to 2 Re^ as found in iHoward et al.l (J2012I ). That rise clearly ends at 2 Re^ The 
occurrence of planets is consistent with constant from 2 Re toward 1 Re- This unexpected plateau in 
planet occurrence at 2 Re suggests distinct planet formation processes for planets above and below 
2 Re' We find that lb.l^\'j% of solar type stars — roughly one in six — has a 1-2 Re planet with 
P = 5-50 days. 
Subject headings: techniques: photometric — planetary systems 



1. INTRODUCTION 

The Kepler Mission has discovered an extraordinary 
sample of more than 2300 planets with radii ranging from 
larg^er than Jupiter to smaller than Earth (Borucki et al.l 
1201 Ibl : iBatalha et all 12012 ). Cleanly measuring and de- 
biasing t his distribution will b e one of Kepler ^s great 
legacies. iHoward et al.l (|2012l ). H12 hereafter, took a 
key step, showing that the planet radius distribution in- 
creases substantially with decreasing planet size down to 
at least 2 Re. While the distribution of planets of all 
periods and radii contains a wealth of information, we 
choose to focus on the smallest planets. Currently, only 
Kepler is able to make quantitative statements about the 
occurrence of planets down to 1 Re- 

The occurrence distributions in HI 2 were based on 
planet candidates detected in the first f our months of 
Kepler photometry (B orucki et al.ll2011b[ ). These planet 
candidates were detected by a sophisticated pipeline de- 
velo ped by the Kepler tediin. Science Operat ions Cen- 
ter (JTwicken et ani2010l : [Jenkins et al.ll2010a[ )Fl Under- 



|epetigura@berkeley.edu | 

^ Since H12. Ba talha et al.l (|2012l ) added many candidates, 
bringing the number of public KOIs to > 2300. In addition. 



standing pipeline completeness, the fraction of planets 
missed by the pipeline as a function of size and pe- 
riod, is a key component to measuring planet occurrence. 
Pipeline completeness can be assessed by injecting mock 
dimmings into photometry and measuring the rate at 
which injected signals are found. The completeness of 
the official Kepler pipeline has yet to be measured in 
this manner. This was the key reason why H12 were 
cautious interpreting planet occurrence under 2 Re- 

In this work, we focus on determining the occurrence 
of small planets. To maximize our sensitivity to small 
planets, we restrict our stellar sample to include only 
the 12,000 stars having the lowest photometric noise in 
the Kepler survey. We comb through quarters 1-12 (Ql- 
Q12) — 3 years of Kepler photometry — with a new al- 
gorithm, TERRA, optimized to detect low signal-to-noise 
transit events. We determine TERRA's sensitivity to plan- 
ets of different periods and radii by injecting synthetic 
transits into Kepler photometry and measuring the re- 
covery rate as a function of planet period and radius. 

We describe our selection of 12,000 low- noise targets 

the Kepler t eam plan et search pipeline has continued to evolve 
(|Smith et a l. 2012 : S tumpe et al. 20121 . 
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in Section [2l We comb their photometry for exoplanet 
transits with TERRA, introduced in Section [3l We re- 
port candidates found with TERRA (Section H]), which 
we combine with our measurement of pipehne complete- 
ness (Section [5j) to produce debiased measurements of 
planet occurrence (Section [6]). We offer some compar- 
isons between TERRA planet candidates and those from 
iBatalha et alJ (|2Q12l ) in Section [3 as well as occurrence 
measured using both catalogs in Section [HI We offer some 
interpretations of the constant occurrence rate for plan- 
ets smaller than 2 Re in Section [9l 

2. THE BEST12K STELLAR SAMPLE 

We restrict our study to the best 12,000 solar type stars 
from the perspective of detecting transits by Earth-size 
planets, hereafter, the "Best 12k" sample. For the small- 
est planets, uncertainty in the occurrence distribution 
stems largely from pipeline incompleteness due to the 
low signal-to-noise ratio (SNR) of an Earth-size transit. 

Our initial sample begins with the 102,835 stars that 
were observed during every quarter from Ql-Q9o From 
this sample, following HI 2, we select 73,757 "solar sub- 
set" stars that are solar-type G and K having Teff = 
4100-6100 K and log^ = 4.0-4.9 (cgs). Teff and log^ 
values are present in the Kepler Input Catalog (KIC; 
iBrown et al.ll201J] ) which is available onlineO Figure [1] 
shows the KIC-based Teff and log^ values as well as 
the solar subset. KIC stellar parameters have large un- 
certainties: aflog.Q ) ~ 0.4 dex and a (Teff) ~ 200 K 
(|Brown et al.ll201lf ). As we will discuss in Section IH we 
determine stellar parameters for the majority of TERRA 
planet candidates spectroscopically. For the remaining 
cases, we use stellar parameters that were determined 
photometrically, but incorporated a main sequence prior 
(IBa talha et al .1120121 ). After refining the stellar parame- 
ters, we find that 10 of the 129 TERRA planet candidates 
fan outside of the Teff = 4100-6100 K and log^ = 4.0- 
4.9 (cgs) solar subset. 

From the 73,757 stars that pass our cuts on log^ and 
Teff, we choose the 12,000 lowest noise stars. Kepler 
target stars have a wide range of noise properties, and 
there are several ways of quantifying photometric noise. 
The Kepler team computes quantities called CDPP3, 
CDPP6, and CDPP12, which are measures of the pho- 
tometric scatter in 3, 6, and 12 hour bins (Jenkins et al. 
l2010a[ ). Since CDPP varies by quarter, we adopt the 
maximum 6-hour CDPP over Q1-Q9 as our nominal 
noise metric. We use the maximum noise level (as op- 
posed to median or mean) because a single quarter of 
noisy photometry can set a high noise floor for planet 
detection. One may circumvent this problem by re- 
moving noisy regions of photometry, which is a planned 
upgrade to TERRA. Figure [2] shows the distribution of 
max(CDPP6) among the 73,757 stars considered for our 
sample. 

In choosing our sample, we wanted to include stars 
amenable to the detection of planets as small as 1 Re- 
We picked the 12,000 quietest stars based on prelimi- 
nary completeness estimates. The noisiest star in the 
Bestl2k sample has max(CDPP6) of 79.2 ppm. We es- 

^ We ran TERRA on Q1-Q12 photometry, but we selected the 
Best 12k sample before Q10-Q12 were available. 
^ http://archive.stsci.edu/Kepler/kic.html 
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Figure 1. Kepler target stars observed every quarter from 
Q1-Q9. The rectangle marks the "solar subset" of stars with 
TefT = 4100-6100 K and log^ = 4.0-4.9 (cgs). 

timated that the ~ 100 ppm transit of an Earth-size 
planet would be detected at SNRcdpp ^ 1.250 Given 
that Q1-Q12 contains roughly 1000 days of photometry, 
we expected to detect a 5-day planet at SNRcdpp ^ 
1.25 X ylOOO/5 ~ 18 (a strong detection) and to detect 

a 50-day planet at SNRcdpp ^ 1-25 x ^"1000/50 - 5.6 (a 
marginal detection) . In our detailed study of complete- 
ness, described in Section [5l we find that TERRA recovers 
most planets down to 1 Re having P = 5-50 days. 

We draw stars from the HI 2 solar subset for two rea- 
sons. First, we may compare our planet occurrence to 
that of H12 without the complication of varying occur- 
rence with different stellar types. We recognize that sub- 
tle differences may exist between the HI 2 and Best 12k 
stellar sample. One such difference is that the Best 12k 
is noise-limited, while the H12 sample is magnitude- 
limited. HI 2 included bright stars with high photomet- 
ric variability, which are presumably young and/or active 
stars. Planet formation efficiency could depend on stellar 
age. Planets may be less common around older stars that 
formed before the metallicity of the Galaxy was enriched 
to current levels. This work assesses planet occurrence 
for a set of stars that are systematically selected to be 
3-10 Gyr old by virtue of their reduced magnetic activity. 

The second reason for adopting the HI 2 solar subset is 
a practical consideration of our completeness study. As 
shown in Section [5l we parameterize pipeline efficiency as 
a function of P and Rp. Because M-dwarfs have smaller 
radii than G-dwarfs, an Earth-size planet dims an M- 
dwarf more substantially and should be easier for TERRA 
to detect. Thus, measuring completeness as a function of 
P, Rp^ and R^ (or perhaps P and Rp/ R^) is appropriate 
when analyzing stars of significantly different sizes. Such 
extensions are beyond the scope of this paper, and we 
consider stars with R^ ~ Rq. 

3. PLANET SEARCH PIPELINE 

Identifying the smallest transiting planets in Kepler 
photometry requires a sophisticated automated pipeline. 
Our pipeline is called "TERRA" and consists of three major 

^ SNRcdpp 5 the expected SNR using the niax(CDPP6) metric, 
is different from the SNR introduced in Section 13.21 SNRcdpp 
is more similar to the SNR computed by the Kepler team, which 
adopts SNRcdpp > 7.1 as their detection threshold. 
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Figure 2. Stellar photometric noise level plotted against Kepler 
magnitude. Noise level is the maximum value of CDPP6 over Ql- 
Q9. Of the 73,757 stars that pass our cuts on Teff and loggf, we 
select the 12,000 most quiet stars. The line shows max(CDPP6) = 
79.2 ppm, corresponding to the noisiest star in the Bestl2k sample, 
well below the median value of 143 ppm. 

components. First, TERRA calibrates photometry in the 
time domain. Then, TERRA combs the cahbrated pho- 
tometry for periodic, box-shaped signals by evaluating 
the signal-to- noise ratio (SNR) over a finely- spaced grid 
in transit period (P), epoch (to) and dur ation (AT). Fi- 
nally, TERRA fits promising signals with a lMandel fc Agoll 
(f2QQ2i) transit model and rejects signals that are not con- 
sistent with an exoplanet transit. We review the calibra- 
tion comp onent in Secti on 13. ![ but refer the reader to 
iPetigura^fc Marcv] (|2Q12l ) for a detailed description. We 
present, for the first time, the grid- sear ch a nd light curve 
fitting components in Sections [X2l and [331 

3.1. Photometric Calibration 

We briefiy review the major time domain components 
of TERRA; for a more comp lete description, please refer to 
iPetigur a fc Marcvl (|2Q12[ ). We begin with Kepler "sim- 
ple aperture long cadence photometry," which we down- 
loaded from the Mikulski Archive for Space Telescopes 
(MAST). This photometry is the total photoelectrons 
accumulated within a predefined target aperture o ver a 
29.4 minute interval (jFraauem fc Thompsonll2Q12[ ). We 
remove thermal settling events manually and cosmic rays 
using a median filter. Next, we remove photometric 
trends longer than 10 days with a high-pass filter. Fi- 
nally, we identify photometric modes shared by a large 
ensemble of stars with using a robust principal compo- 
nents analysis. The optimum linear combination of the 
four most significant modes is removed from each light 
curve individually. 

3.2. Grid-based Transit Search 

We then search for periodic, box-shaped signals in 
ensemble-calibrated photometry. Such a search involves 
evaluating the SNR over a finely sampled grid in period 



(P), epoch (to), and duration (AT), i.e. 

SNR = SNR(P,to,AT). (1) 

O ur approach is simila r to the widely-used BLS algorithm 
of lKovacs et al.l (|2QQ2l) as well as to the TP S component of 
the Kepler pipeline (| Jenkins et al.ll2QlQb ). BLS, TPS, and 
TERRA are all variants of a "matched filter" (North 1943). 
The way in which such an algorithm searches through 
P, to, and AT is up to the programmer. We choose to 
search first through AT (outer loop), then P, and finally, 
to (inner loop). 

For computational simplicity, we consider transit du- 
rations that are integer numbers of long cadence mea- 
surements. Since we search for transits with P = 5-50 
days, we try AT = [3, 5, 7, 10, 14, 18] long cadence mea- 
surements, which span the range of expected transit du- 
rations, 1.5 to 8.8 hours, for G and K dwarf stars. 

After choosing AT, we compute the mean depth, 
5F{ti)^ of a putative transit with duration = AT cen- 
tered at ti for each cadence. 5F is computed via 



6F{U) = Y.P{U-i)^, 



(2) 



where F(ti) is the median- normalized stellar fiux at time 
ti and Gj is the jth element of the following kernel 
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(3) 



(4) 



We search over a finely sampled grid of trial peri- 
ods from 5-50 days and epochs ranging from tgtart to 
4tart + Pi where tgtart is the time of the first photometric 
observation. For a given (P^to^AT) there are Nt puta- 
tive transits with depths 6Fi, for i = 0, 1, ... , Nt — 1. 
For each (P,to,AT) triple, we compute SNR from 



SNR: 



-mean((5Pi), 



(5) 



where a is a robust estimate (median absolute deviation) 
of the noise in bins of length AT. 

For computational efficiency, we empl oy the "Fast Fold- 
ing Algorithm" (FFA) of St aelinl (|1969f ) as implemented 
in Petigura fc Marcy (2013; in prep.). Let Pcad,o be 
a trial period that is an integer number of long ca- 
dence measurements, e.g. Pcad,o = 1000 implies P = 
1000 X 29.4 min = 20.43 days.' Let A^cad = 51413 be 
the length of the Q1-Q12 time series measured in long 
cadences. Leveraging the FFA, we compute SNR at the 
following periods: 



cad,i 



cad.O 



M-1' 



0,1,...,M-1 (6) 



where M = A'cad/^cad,o rounded up to the nearest power 
of two. In our search from 5-50 days, Pcad,o ranges from 
245-2445, and we evaluate SNR at ~ 10^ different peri- 
ods. At each Pcad i we evaluate SNR for Pcad o different 
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starting epochs. All told, for each star, we evaluate SNR 
at ~ 10^ different combinations of P, to, and AT. 

Due to runtime and memory constraints, we store only 
one SNR value for each of the trial periods. TERRA stores 
the maximum SNR at that period for all AT and to. 
We refer to this one-dimensional distribution of SNR as 
the "SNR periodogram," and we show the KIC-3120904 
SNR periodogram in Figure [3] as an example. Because 
we search over many AT and to at each trial period, 
fluctuations often give rise SNR ~ 8 events and set the 
detectability floor in the SNR pe riodogram. Fo r KIC - 
3120904, a star not listed in the iBatalha et all (|2012r ) 
planet catalog, we see a SNR peak of 16.6, which rises 
clearly above stochastic background. 

If the maximum SNR in the SNR periodogram ex- 
ceeds 12, we pass that particular (P,to,AT) on to the 
"data validation" (DV) step, described in the following 
section, for additional vetting. We chose 12 as our SNR 
threshold by trial and error. Note that the median ab- 
solute deviation of many samples drawn from a Gaus- 
sian distribution is 0.67 times the standard deviation, 
i.e. CTMAD = 0.67crsTD. Therefore, TERRA SNR = 12 
corresponds roughly to SNR = 8 in a BLS or TPS search. 

Since TERRA only passes the (P,to,AT) triple with the 
highest SNR on to DV, TERRA does not detect additional 
planets with lower SNR due to either smaller size or 
longer orbital period. As an example of TERRA's insen- 
sitivity to small candidates in multi-candidate systems, 
we show the TERRA SNR periodogram for KIC-5094751 
in Figure [H Batalha et al. (2012) lists two candidates 
belonging to KIC-5094751: KOI-123.01 and KOI-123.02 
with P = 6.48 and 21.22 days, respectively. Although 
the SNR periodogram shows two sets of peaks coming 
from two distinct candidates, TERRA only identifies the 
first peak. Automated identification of multi-candidate 
systems is a planned upgrade for TERRA. Another caveat 
is that TERRA assumes strict periodicity and struggles to 
detect low SNR transits with significant transit timing 
variations, i.e. variations longer than the transit dura- 
tion. 

3.3. Data Validation 

If the SNR periodogram has a maximum SNR peak > 
12, we flag the corresponding (P,to,AT) for additional 
vetting. Following the language of the official Kepler 
pipeline, we refer to these triples as "threshold crossing 
events" (TCEs), since they have high photometric SNR, 
but are not necessarily consistent with an exoplanet tran- 
sit. TERRA vets the TCEs in a step called "data valida- 
tion," again following the nomenclature of the official Ke- 
pler pipeline. Data validation (DV), as implement ed in 
the official Kepler pipeline, is described in Jenkins et al.l 
(j2010ar ). We emphasize that TERRA DV does not depend 
on the DV component of the Kepler team pipeline. 

We show the distribution of maximum SNR for each 
Bestl2k star in Figured Among the Bestl2k stars, 738 
have a maximum SNR peak exceeding 12. Adopting SNR 
= 12 as our threshold balances two competing needs: the 
desire to recover small planets (low SNR) and the desire 
to remove as many non-transit events as possible before 
DV (high SNR). As discussed below, only 129 out of ah 
738 events with SNR > 12 are consistent with an exo- 
planet transit, with noise being responsible for the re- 
maining 609. As shown in Figure [5l that number grows 



rapidly as we lower the SNR threshold. For example, 
the number of TCEs grows to 3055 with a SNR thresh- 
old of 10, dramatically increasing the burden on the DV 
component. 

A substantial number (347) of TCEs are due to har- 
monics or subharmonics of TCEs outside of the P = 5-50 
day range and are discarded. In order to pass DV, a TCE 
must also pass a suite of four diagnostic metrics. The 
metrics are designed to test whether a light curve is con- 
sistent with an exoplanet transit. We describe the four 
metrics in Table [1] along with the criteria the TCE must 
satisfy in order pass DV. The metrics and cuts were de- 
termined by trial and error. We recognize that the TERRA 
DV metrics and cuts are not optimal and discard a small 
number of compelling exoplanet candidates, as discussed 
in Section [731 However, since we measure TERRA's com- 
pleteness by injection and recovery of synthetic transits, 
the sub-optimal nature of our metrics and cuts is incor- 
porated into our completeness corrections. 

Our suite of automated cuts removes all but 145 TCEs. 
We perform a final round of manual vetting and re- 
move 16 additional TCEs, leaving 129 planet candidates. 
Most TCEs that we remove manually come from stars 
with highly no n- stationary photometric noise properties. 
Some stars have small regions of photometry that exceed 
typical noise levels by a factor of 3. We show the SNR 
periodogram for one such star, KIC-7592977, in Figure[6l 
Our definition of SNR (Equation [5j) incorporates a sin- 
gle measure of photometric scatter based on the median 
absolute deviation, which is insensitive to short bursts 
of high photometric variability. In such stars, fluctua- 
tions readily produce SNR ~ 12 events and raise the 
detectability floor to SNR - 12, up from SNR - 8 in 
most stars. We also visually inspect phase-folded light 
curves for coherent out-of-transit variability, not caught 
by our automated cuts, and for evidence of a secondary 
eclipse. 

4. SMALL PLANETS FOUND BY TERRA 

Out of the 12,000 stars in the Bestl2k sample, TERRA 
detected 129 planet candidates achieving SNR > 12 that 
passed our suite of DV cuts as well as visual inspection. 
Table [2] lists the 129 planet candidates. We derive planet 
radii using Rp/R^ (from Mandel-Agol model fits) and R^ 
from spectroscopy (when available) or broadband pho- 
tometry. 

We obtained spectra for 100 of the 129 stars using 
HIRES (Vogt et al. 1994) at the Keck I telescope with 
standard configura tion of the California Planet Survey 
(jMarcv et al.l 120081 ). These spectra have resolution of ~ 
50,000, at a signal-to-noise of 45 per pixel at 5500 A. 
We determine stellar parameters using a routine called 
SpecMatch (Howard et al. 2013, in prep). In brief, 
SpecMatch compares a stellar spectrum to a library of 
- 800 spectra with Teff = 3500-7500 K and log^ = 2.0- 
5.0 (determined from LTE spectral modeling). Once the 
target spectrum and library spectrum are placed on the 
same wavelength scale, we compute x^^ the sum of the 
squares of the pixel- by-pixel differences in normalized in- 
tensity. The weighted mean of the ten spectra with the 
lowest x^ values is taken as the final value for the effec- 
tive temperature, stellar surface gravity, and metallicity. 
We estimate SpecMatch-derived stellar radii are uncer- 
tain to 10% RMS, based on tests of stars having known 
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Figure 3. SNR periodogram of KIC-3120904 photometry. We evaluate SNR over a finely-spaced, three-dimensional grid of P, to, and 
AT. We store the maximum SNR for each trial period, resulting in a one-dimensional distribution of SNR. A planet candidate (not in 
IBatalha et al.ll20T^ produces a SNR peak of 16.6 at P = 42.9 days, which rises clearly above the detection fioor of SNR ~ 8. 
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Figure 4. SNR periodogram of KIC-5094751 photometry, demonstrating TERRA's insensitivity to lower SNR candidates in multi-candidate 
systems. Batalha et al. (2012) lists two planets belonging to KIC-5094751, KOI-123.01 and KOI-123.02 with P = 6.48 and 21.22 days, 
respectively. TERRA detected KOI-123.01 with a period of 6.48 days (highest SNR peak). Sub-harmonics belonging to KOI-123.01 are 
visible at [2, 3, . . .] x P = [13.0, 19.4, . . .] days. A second set of SNR peaks due to KOI-123.02 (P = 21.2 days) is visible at [0.5, 2, . . .] x P = 
[10.6,42.4, . . .] days. Had we removed the transit due to KOI-123.01, KOI-123.02 would be easily detectible due its high SNR of ~ 80. 
TERRA does not yet include multi-candidate logic and is thus blind to lower SNR candidates in multi-candidate systems. 



In order to pass the "data validation' 



Table 1 

(DV) stage, a "threshold crossing event" 
of cuts. 



(TCE) must pass the following suite 



description 



value 



s2n_out_on_in 



med_on_mean 



Compelling transits have fiat out-of-transit light curves. For a TCE with (P,to,AT), 
we remove the transit region from the light curve and evaluate the SNR of all other 
{P' ,t'Q,AT') triples where P = P' and AT = AT^ s2n_out_on_in is the ratio of the 
two highest SNR events. 

Since the our definition of SNR (Equation [5]) depends on the arithmetic mean of indi- 
vidual transit depths, outliers occasionally produce high SNR TCEs. For each TCE, we 
compute a robust SNR, 



< 0.7 



> 1.0 



medSNR = 



- median ((^P^). 



med_on_mean is medSNR divided by SNR as defined in Equation [5] 

We compute the circular autocorrelation of the phase-folded light curve, autor is the 

ratio of the highest autocorrelation peak (at lag) to the second highest peak and is 

sensitive to out-of-transit variability. 

We fit the phase-folded light curve with a IMandel &: Agoll (|2002l ) model, taur is the 
ratio of the best fit transit duration to the maximum duration given the KIC stellar 
parameters and assuming a circular orbit. 



> 1.6 



< 2.0 



radii from high resolution spectroscopy and asteroseis- 
mology. 

For 27 stars where spectra are not available, we 
adopt the pho tomet rically- derived stellar parameters of 
[Batalha et aLi ([20121). These parameters are taken from 
the KIC (jBrown et all 1201 If ), but then modified so that 
th ey lie on the Y onsei-Yale stellar evolution models 
of iDemarque et aTl (J2004). The resulting stellar radii 
have uncertainties of 35% (rms), but can be incor- 



rect by a factor of 2 or more. As an extreme ex- 
ample, the interpretati ons of the thr ee planets in the 
KOI-961 system (Muirh ead et al.|[2012l ) changed dramat- 
ically when HIRES spectra showed the star to be an 
M5 dwarf (0.2 i?© as opposed to 0.6 Rq listed in the 
KIC). We could not obtain spectra for two stars, KIC- 
7345248 and KIC-84 29668, which were not present in 
IBatalha et al.l (|2012[ ). We determine stellar parameters 
for these stars by fitting the KIC photometry to Yonsei- 
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Figure 5. Distribution of the highest SNR peak for each star 
in the Bestl2k sample. We show SNR = 5-20 to highhght the 
distribution of low SNR events. The 738 stars with SNR > 12 are 
labeled "threshold crossing events" (TCEs) and are subjected to 
additional scrutiny in the "data validation" component of TERRA. 



Yale stellar models. We adopt 35% fractional errors on 
photometrically-derived stellar radii. 

On ce we determine P and to, we fit a'Mand el &: Agoll 
(|2QQ2 ) model to the phase-folded photometry. Such a 
model has three free parameters: Rp/R^, the planet to 
stellar radius ratio; r, the time for the planet to travel 
a distance R^ during transit; and 6, the impact param- 
eter. In this work, Rp/R^ is the parameter of interest. 
However, b and Rp/R^ are covariant, i.e. a transit with 
b approaching unity only traverses the limb of the star, 
and thus produces a shallower transit depth. In order 
to account for this covariance, best fit parameters were 
computed via Markov Chain Monte Carlo. We find that 
the fractional uncertainty on Rp/R^^ ^i^ ^/r "^ ^^^ ^^ ^^ 
high as 10%, but is generally less than 5%. Therefore, 
the error on Rp due to covariance with b is secondary to 
the uncertainty on R^. 

We show the distribution of TERRA candidates in Fig- 
ure [7] over the two-dimensional domain of planet radius 
and orbital period. Our 129 candidates range in size from 
6.83 Re to 0.48 Re (smaller than Mars). The median 
TERRA candidate size is 1.58 Re- In Figure [HI we show 
the substantial overlap between the TERRA planet sample 
and those produced by the Ke pler te a m. TE RRA recovers 
82 candidates listed in Batalh a et al.l ()2012f ) . We discuss 
the significant overlap between the two works in detail 
in Section As of August 8th, 2012, 10 of our TERRA 
candidates were listed as false positives in an internal 
database of Kepler planet candidates maintained by Ja- 
son Rowe (Jason Rowe, 2012, private communication) 
and are shown as blue crosses in Figure [51 We do not in- 
clude these 10 candidates in our subsequent calculation 
of occurrence. Table [2] lists the KIC identifier, best fit 
transit parameters, stellar parameters, planet radius, and 
Kepler team false positive designation of all 129 candi- 
dates revealed by the TERRA algorithm. The best fit tran- 
sit parameters include orbital period, P; time of transit 
center, to; planet to star radius ratio, Rp/R^; time for 
planet to cross R^ during transit, r; and impact param- 
eter, b. We list the following stellar properties: effective 
temperature, Tgff; surface gravity, log^; and stellar ra- 



dius, Ri,. 

5. COMPLETENESS OF PLANET CATALOG 

When measuring the distribution of planets as a func- 
tion of P and Rp, understanding the number of missed 
planets is as important as finding planets themselves. 
HI 2 accounted for completeness in a rough sense based 
on signal-to- noise considerations. For each star in their 
sample, they estimated the SNR over a range of P and 
Rp using CDPP as an estimate of the photometric noise 
on transit-length timescales. H12 chose to accept only 
planets with SNR > 10 in a single quarter of photome- 
try for stars brighter than Kp = 15. This metric used 
CDPP and was a reasonable pass on the data, particu- 
larly when the pipeline completeness was unknown. De- 
termining expected SNR from CDPP does not incorpo- 
rate the real noise characteristics of the photometry, but 
instead approximates noise on transit timescales as sta- 
tionary (CDDP assumed to be constant over a quarter) 
and Gaussian distributed. Moreover, identifying small 
transiting planets with transit depths comparable to the 
noise requires a complex, multistage pipeline. Even if 
the integrated SNR is above some nominal threshold, the 
possibility of missed planets remains a concern. 

We characterize the completeness of our pipeline by 
performing an extensive suite of injection and recovery 
experiments. We inject mock transits into raw photom- 
etry, run this photometry though the same pipeline used 
to detect planets, and measure the recovery rate. This 
simple, albeit brute force, technique captures the idiosyn- 
crasies of the TERRA pipeline that are missed by simple 
signal-to-noise considerations. 

We perform 10,000 injection and recovery experiments 
using the following steps: 

1. We select a star randomly from the Best 12k sam- 
ple. 

2. We draw (P^Rp) randomly from log-uniform dis- 
tributions over 5-50 days and 0.5-16.0 Re- 

3. We draw impact parameter and orbital phase ran- 
domly from uniform distributions ranging from 
to 1. 

4. We generate a lMandel fc Agoll (I2002[ ) model. 

5. We inject it into the "simple aperture photometry" 
of the selected star. 

We then run the calibration, grid-based search, and data 
validation components of TERRA (Sections 13. If 13. 2[ and 
13. 3p on this photometry and calculate the planet recovery 
rate. We do not, however, perform the visual inspection 
described in Section[331 An injected transit is considered 
recovered if the following two criteria are met: (1) The 
highest SNR peak passes all DV cuts and (2) the output 
period and epoch are consistent with the injected period 
and epoch to within 0.01 and 0.1 days, respectively. 

Figure [9] shows the distribution of recovered simula- 
tions as a function of period and radius. Nearly all simu- 
lated planets with Rp > 1.4 Re are recovered, compared 
to almost none with Rp < 0.7 Re- Pipeline completeness 
is determined in small bins in (P,i?p )-space by dividing 
the number of successfully recovered transits by the total 
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exhibited short bursts of high photometric scatter, which raised the noise floor to SNR ~ 12, up from SNR ~ 8 as in most stars. 
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Figure 7. Periods and radii of 129 planet candidates detected by TERRA. Errors on Rp are computed via 



Rp 



^ ^/p *^ j , where Rp/Ri, is the radius ratio. The error in Rp stems largely from the uncertainty in stellar radii. 



We adopt 



R* 



10% for the 100 stars with spectroscopically determined R-^ and 



r(R.) _ 



35% for the remaining stars with R^, deter- 



mined from photometry. Using MCMC, we flnd the uncertainty in Rp / R^, is generally < 5% and thus a minor component of the overall 
error budget. 



number of injected transits in a bin- by-bin basis. This 
ratio is TERRA's recovery rate of putative planets within 
the Bestl2k sample. Thus, our quoted completeness esti- 
mates only pertain to the low photometric noise Best 12k 
sample. Had we selected an even more rarified sample, 
e.g. the "BestGk," the region of high completeness would 
extend down toward smaller planets. 

6. OCCURRENCE OF SMALL PLANETS 

Following H12, we define planet occurrence, /, as the 
fraction of a defined population of stars having planets 
within a domain of planet radius and period, including 
all orbital inclinations. TERRA, however, is only sensitive 
to one candidate (highest SNR) per system, so we report 
occurrence as the fraction of stars with one or more plan- 
ets with P = 5-50 days. Our occurrence measurements 
apply to the Bestl2k sample of low-noise, solar-type stars 
described in Section O 



In computing planet occurrence in the Best 12k sam- 
ple, we follow the prescription in H12 with minor modi- 
fications. Notably, we have accurate measures of detec- 
tion completeness described in the previous section. In 
contrast, H12 estimated completeness based on the pre- 
sumed signal-to- noise of the transit signal, suffering both 
from approximate characterization of photometric noise 
using CDPP and from poor knowledge of the efficiency 
of the planet-finding algorithm for all periods and sizes. 

For each P-Rp bin, we count the number of planet 
candidates, npi^ceii- Each planet that transits represents 
many that do not transit given the orientation of their 
orbital planes with respect to Kepler^ s line of sight. As- 
suming random orbital alignment, each observed planet 
represents a/i?^ total planets when non-transiting ge- 
ometries are considered. For each cell, we compute the 
number of augmented planets, npi^aug,ceii = ^iCa/ R-k,i^ 
which accounts for planets with non-transiting geome- 





8.0 




5.7 




4.0 




2.8 


'-0 

2 


2.0 
1.4 


LU 


1.0 


Q_ 


0.7 



E. Petigura & G. Marcy 



0.5 



■ X 



■ ■ 



TERRA candidates 



82 in Bataiha et al. (2012) ' 

10 listed as false positives ^ 

in Kep/er database (Aug, 2012) 



18 additional candidates 

in Kep/er database (Aug, 2012) 

19 unlisted candidates 



— I — 
10 



23 



50 



Period [days] 
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Figure 9. Results from the injection and recovery of 10,000 synthetic transit signals into actual photometry of randomly selected stars 
from our Bestl2k stellar sample. Each point represents the planet radius and orbital period of a mock transiting planet. The blue points 
represent signals that passed the DV post-analysis and where TERRA recovers the correct period and epoch. Signals that did not pass DV 
and/or were not successfully recovered, are shown as red points. Pipeline completeness is simply the number of blue points divided by the 
total number of points in each bin. The figure shows that for planet sizes above 1.0 Re, our pipeline discovers over 50% of the injected 
planets, and presumably accomplishes a similar success rate for actual transiting planets. The completeness for planets larger than 1 Re 
is thus high enough to compute planet occurrence for such small planets, with only moderate completeness corrections needed (less than a 
factor of 2). Note that we are measuring the recovery rate of putative planets in the Best 12k sample with TERRA. Had we selected a lower 
noise stellar sample, for example the "BestGk," the region of high completeness would extend to even small radii. 



tries. We then use Kepler's 3^^ law together with P and 
M^ to compute a/i?^ assuming a circular orbit 

To compute occurrence, we divide the number of stars 
with planets in a particular cell by the number of stars 
amenable to the detection of a planet in a given cell, 
^•,amen- This uumbcr is just N^ = 12,000 times the 
completeness, computed in our Monte Carlo study. The 
debiased fraction of stars with planets per P-Rp bin, 

/cell, is given by /cell = npi,aug,cell/n^,amen. We shoW /cell 

on the P-Rp plane in Figure [TU] as a color scale. We 
also compute (Pfce\\/d\ogP/d\ogRp^ i.e. planet occur- 
rence divided by the logarithmic area of each cell, which 
is a measure of occurrence which does not depend on 
bin size. We annotate each P-Rp bin of Figure [TOl with 
the corresponding value of npi,ceii, npi,aug,ceii, /cell, and 
d^fceii/dlogP/dlogRp. 

Due to the small number of planets in each cell, er- 
rors due to counting statistics alone are significant. We 
compute Poisson errors on npi^ceii for each cell. Errors on 
^pi,aug,ceii, /cell, and d'^ fcew/ dlog P/ dlog Rp include only 
the Poisson errors from npi^ceii- There is also shot noise 
associated with the Monte Carlo completeness correction 
due to the finite number of simulated planets in each P- 
Rp cell, but such errors are small compared to errors on 
^pi,ceii- The orbital alignment correction, a/R^, is also 
uncertain due to imperfect knowledge of stellar radii and 
orbital separations. We do not include such errors in our 
occurrence estimates. 

Of particular interest is the distribution of planet oc- 
currence with Rp for all periods. We marginalize over 
P by summing occurrence over all period bins from 5 to 
50 days. The distribution of radii shown in Figure [TTI 
shows a rapid rise in occurrence from 8.0 to 2.8 Re- 
HI 2 also observed a rising occurrence of planets down 
to 2.0 Re^ which they modeled as a power law. Planet 
occurrence is consistent with a flat distribution from 2. 8 
to 1.0 Re, ruling out a continuation of a power law in- 
crease in occurrence for planets smaller than 2.0 Re- 
We find lb.lt_\j% of Sun-like stars harbor a 1.0-2.0 Re 
planet with P = 5-50 days. Including larger planets, we 
find that 24.8134% of stars harbor a planet larger than 
Earth with P = 5-50 days. Occurrence values assum- 
ing a 100% efficient pipeline are shown as gray bars in 
Figure [TTI The red bars show the magnitude of our com- 
pleteness correction. Even though TERRA detects many 
planets smaller than 1.0 Re^ we do not report occurrence 
for planets smaller than Earth since pipeline complete- 
ness drops abruptly below 50%. 

We show planet occurrence as a function of orbital pe- 
riod in Figure [121 In computing this second marginal 
distribution, we include radii larger than 1 Re so that 
corrections due to incompleteness are small. Again, as 
in Figure [m gray bars represent uncorrected occurrence 
values while red bars show our correction to account for 
planets that TERRA missed. Planet occurrence rises as 
orbital period increases from 5.0 to 10.8 days. Above 
10.8 days, planet occurrence is nearly constant per loga- 
rithmic period bin with a slight indication of a continued 
rise. This leveling off of the distribution was noted by 
H12, who considered Rp > 2.0 Re- We fit the distribu- 

^ H12 determined a/i?^ directly from light curve fits, but found 
little change when computing occurrence from a/i?^ using Kepler's 
third law. 



tion of orbital periods for Rp > 1.0 Re with two power 
laws of the form 



df 



dlogP 



kpP'', 



(7) 



where a and kp are free parameters. We find best fit 
values of kp = 0.185i^;^^^, a = 0.16 ± 0.07 for P = 5- 

10.8 days and kp = S.Atol x 10-^ a = 1.35 ± 0.05 for 
P = 10.8-50 days. We note that kp and a are strongly 
covariant. Extrapolating the latter fit speculatively to 
P > 50 days, we find 41.7t5;9% ^^ Sun-like stars host a 
planet 1 Re or larger with P = 50-500 days. 

7. COMPARISON OF TERRA AND BATALHA ET AL. (2012) 
PLANET CATALOGS 

Her e, we compare our candidates to those of 
Batalha et al.l (2012). Candidates were deemed in com- 
mon if their periods agree t o within 0.01 days. W e list 
the union of the TERRA and iBatalha et al.l (|2012f ) cata- 
logs in Table HI Eighty-two candidates appear in both 
catalogs (Section l7.ip . 47 appear in this work only (Sec- 
tion 17. 2p . and 33 appear in Batalha et al. (2012) only 
(Section [7]3j). We discuss the significant overlap between 
the two catalogs and explain why some candidates were 
detected by one pipeline but not the other. 

7.1. Candidates in Common 

Eighty-two of our candidates appear in the 

'Batalh a et al.l (1201 2*) catalog. We show these can- 
didates in P-Rp space in Figure [8] as grey points. 
TERRA detected no new candidates with Rp > 2 Re- 
This agreement in detected planets having Rp > 2 Re 
demonstrates high completeness for such planets in both 
pipelines for this sample of quiet stars. This is not very 
surprising since candidates with Rp > 2 Re have high 
SNR, e.g. min, median, and max SNR = 19.3, 71.5, and 
435 respectively. 

Radii for the 82 planets in c o mmon were fairly con- 
sistent between IBatalha et al.l ()2012f ) and this work. 
The two exceptions were KIC-8242434 and KIC-8631504. 
Using SpecMatch, we find stellar radii of 0.68 and 
0.72 R(r), respectively, down from 1.86 and 1.80 Rq 
in Batalha et al] (|2012i! ). The revised planet radii are 
smaller by over a factor of two. Radii for the other plan- 
ets in common were consistent to ~ 20%. 



7.2. TERRA Candidates Not in Watalha et al.\ 11201^ ) 
Catalog 

TERRA reveal ed 47 p l anet c andidates that did not ap- 
pear in Batalh aet al.l (|2012l ). Such candidates are col- 
ored blue and red in Figure [51 Many of these new de- 
tections likely stem from the fact that we use twice the 
photometry that was available to Batalha et al. (2012). 
To get a sense of how additional photometr y improves the 
planet yield of the Kepler pipeline bevond lBatalha et al.l 
([20121), we compared the TERRA candidates to the Kepler 
team KOI list dated August 8, 2012 (Jason Rowe, private 
communication). The 28 candidates in common between 
the August 8, 2012 Kepler team sample and this work 
are colored blue in Figure [51 Of these 28 candidates, 
10 are listed as false positives and denoted as crosses in 
Figure [51 
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Figure 10. Planet occurrence as a function of orbital period and planet radius for P = 5-50 days and Rp = 0.5-8 Re- TERRA planet 
candidates are shown as red points. Cell occurrence, /cell 5 is given by the color scale. We quote the following information for each cell: 
Top left-number of planets (number of augmented planets); lower left-completeness; top right-fractional planet occurrence, /cell; bottom 
right-normalized planet occurrence, d? fce\\/d\ogP/d\ogRp. We do not color cells where the completeness is less than 50% (i.e. the 
completeness correction is larger than a factor of 2). 

We announce 3 7 ne w planet candidates with respect 
to'B atalha et alJ (^2012) that were not hsted as false pos- 
itives in the Kepler team sample. These 37 candidates, 
all with Rp ^ 2 Re , are a subset of those listed in Ta- 
ble [21 As a convenience, we show this subset in Table [3l 
We remind the reader that all photometry used in this 
work is publicly available. We hope that interested read- 
ers will fold the photometry on the ephemeris in Table [2] 
and access critically whether a planet interpretation is 
correct. As a quick reference, we have included plots 
of the transits of the 37 new candidates from Table [3] 
in the appendix (Figures [TTl and [T8|) . We do not claim 
that our additional candidates bring pipeline complete- 
ness to unity for planets with Rp ^ 2 Re- As shown 
in Section [5l our planet sample suffers from significant 
incompleteness in the same P-Rp space where most of 
the new candidates emerged. 



7.3. \Batalha ef 



11201^ ) Candidates Not in TERRA 
catalog 



The re are 33 planet candidates in the iBatalha et al.l 
(j2012f ) catalog from Bestl2k stars that TERRA missed. Of 
these, 28 are multi-candidate systems where one compo- 
nent was identified by TERRA. TERRA is currently insen- 
sitiv e to multiple planet systems (as des cribed in Sec- 
tion 13. 2p . TERRA missed the remaining 5 IBatalha et al.l 
(|2012l ) candidates for the following reasons: 

• 2581.01 : A bug in the pipeline prevented success- 
ful photometric calibration (Section l3.ip . This bug 
affected 19 out of 12,000 stars in the Best 12k sam- 
ple. 

• 70.01, 111.01, 119.01 : Failed one of the automated 
DV cuts (taur, med_on_mean, and taur, respec- 
tively). We examined these three light curves in 
the fashion described in Section 13. 3[ and we deter- 
mined these light curves were consistent with an 
exoplanet transit. The fact that DV is discarding 
compelling transit signals decreases TERRA's overall 
completeness. Computing DV metrics and choos- 
ing the optimum cuts is an art. There is room for 
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Figure 11. Distribution of planet occurrence for Rp ranging from 
1.0 to 8.0 Re- We quote the fraction of Sun-like stars harboring a 
planet with P = 5-50 days for each Rp bin. We observe a rapid 
rise in planet occurrence from 8.0 down to 2.8 Re, as seen in H12. 
Below 2.8 Re, the occurrence distribution is consistent with flat. 
This result rules out a power law increase in planet occurrence to- 
ward smaller radii. Adding up the two smallest radius bins, we 
flnd 15.1_2'7% of Sun-like stars harbor a 1.0-2.0 Re planet within 
~ 0.25 AU. To compute occurrence as a function of Rp, we simply 
sum occurrence rates for all period bins shown in Figure 1101 Er- 
rors due to counting statistics are computed by adding errors from 
each of the three period bins in quadrature. The gray portion of 
the histogram shows occurrence values before correcting for missed 
planets due to pipeline incompleteness. Our correction to account 
for missed planets is shown in red, and is determined by the injec- 
tion and recovery of synthetic transits described in Section [5] We 
do not show occurrence values where the completeness is < 50%. 




10.8 15.8 23.2 
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Figure 12. Distribution of planet occurrence for different orbital 
periods ranging from 5 to 50 days. We quote the fraction of Sun- 
like stars with a planet Earth-size or larger as a function of orbital 
period. We observe a gradual rise in occurrence from 5.0 to 10.8 Re 
followed by a leveling off for longer orbital periods. H12 observed 
a similar leveling off in their analysis which included planets larger 
than 2 Re- We flt the domains above and below 10.8 Re separately 
with power laws, df/d\ogP = kpP^. We flnd best flt values of 



kp = 0.185; 



+0.043 



-0.035' 

a - 



: 0.16 ± 0.07 for P = 5-10.8 days and kp 



8.4t°-|xl0-3, a = 1.35±0.05 for P = 10.8-50 days. Speculatively, 
we extrapolate the latter power law flt another decade in period 
and estimate 41.7^5g% of Sun-like stars harbor a planet Earth- 
size or larger with P = 50-500 days. As in Figure 1111 the gray 
portion of the histogram shows uncorrected occurrence while the 
red region shows our correction for pipeline incompleteness. Note 
that the number of detected planets decreases as P increases from 
10.8 to 50 days, while occurrence remains nearly constant. At 
longer periods, the geometric transit probability is lower, and each 
detected planet counts more toward df/d log P. 



improvement here. 

• K0I-1151.Q1 : Perio d misidentifi ed in 

iBatalha et all (|2012r ). In iBatalha et all (|2012r ) 
KOI-1151.01 is listed with with a P = 5.22 days. 
TERRA found a candidate with P = 10.43 days. 
Figure [13] shows phase- folded photometry with the 
TERRA ephemeris. A period of 5.22 days would 
imply dimmings in regions where the light curve is 
fiat. 



We plot the 33 total candidates listed in IBatalha et al.l 
(j2012f ). but not found by TERRA in Figure [Ml We high- 
light the 5 missed candidates that cannot be explained 
by the fact that they are a lower SNR candidate in a 
multi-candidate system. TERRA is blind to planets in sys- 
tems with another planet with higher SNR. Figure [Ml 
shows that most of these missed planets occur at Rp < 
1.4 Re. 
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Figure 13. Phase-folded photometry of KIC-8280511 folded on 
the correct 10.43 day period found by TERRA. KOI-1151.01 is listed 
with P = 5.22 days in Batalha et al. (2012). If the transit was 
truly on the 5.22 day period, we should see a transit of equal depth 
180 de grees out of phase. KOI-1151.01 is listed in .Batalha et al.1 
(|2012l ) with half its true period. 



8. OCCURRENCE WITH PLANET MULTIPLICITY 
INCLUDED 

While the TERRA planet occurrence measurement ben- 
efits from well-characterized completeness, it does not 
include the contribution of multis t o ov erall pla net oc- 
currence. As discussed in Sections 13.21 and 17. 3[ TERRA 
only detects the highest SNR candidate for a given star. 
Here, we present pla net occurrence including multis from 
IBatalha et al.l (|2Q12l ). Thus, the occurrence within a bin, 
/, in this section should be interpreted as the average 
number of planets per s tar with P = 5-50 days. The 
additional planets from Batalha et al.l (|2Q12l ) raise the 
occurrence values somewhat over those of the previous 
section. However, the rise and plateau structure remains 
the same. 

We comput e /cell from the 32 candidates present in 
IBatalha et aH (|2012l but not found by TERRA (misla- 
beled KOI-1151.01 was not included). For clarity, we re- 
fer to this separate occurrence calculation as /cell, Batalha- 
Because the completeness of the Kepler pipeline is un- 
known, we apply no completeness correction. This as- 
sumption of 100% completeness is certainly an overesti- 
mate, but we believe that the sensitivity of the Kepler 
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Figure 14. P and Rp for the 33 candidates present in I Bataiha et al.l (|2012 ) but not found by TERRA. The small symbols show the 
candidates in mulit-planet systems. TERRA is blind to such candidates. The 5 larger symbols show the other failure modes of TERRA: 2581.01 
failed due to a pipeline bug; 70.01, 111.01, and 119.01 did not pass DV; and TERRA missed KOI-1151.01 because it is listed in lBatalha et al.l 
(f2012. ') with the incorrect period. Most of the missed planets have Rp < 1.4 Re- 



pipeline to multis is nearly complete for Rp > 1.4 Re- 
TERRA has > 80% completeness for Rp > 1.4 Re because 
planets in that size range with P = 5-50 days around 
Best 12k stars have high SNR. The Kepler pipeline should 
also be detecting these high SNR candidates. Also, once 
a KOI is found, the Kepler team reprocesses the light 
curve for additional transits (Jason Rowe, private com- 
munication). Due to this additional scrutiny, we believe 
that the Kepler completeness for multis is higher than 
for singles, all else being equal. 

We then add /ceii,Bataiha to /ceii computed in the previ- 
ous secti on. We show oc currence computed using TERRA 
and Bataiha et al.l (|2012l ) planets as a function of P and 
Rp in Figure [15] and as a function of on ly Rp in Fig- 
ure [TBI The 32 additional planets from iBatalha et al.l 
()2012l ) do not change the overall shape of the occurrence 
distribution: rising from 4.0 to 2.8 Re and consistent 
with flat from 2.8 down to 1.0 Re- 

H12 fit occurrence for Rp > 2 Re with a power law. 



df 



d log Rp 



kpR"^, 



(8) 



finding a = -1.92 ± 0.11 and kp = 2.9^0.1 (Section 
3.1 of H12). As a point of comparison, we plot the 
H12 power law over our combined occurrence distribu- 
tion in Figure [ini The fit agrees qualitatively for Rp > 
2 Re, but not within errors. We expect the H12 fit to be 
^ 25% higher than our occurrence measurements since 
HI 2 included planets with P < 50 days (not P = 5- 
50 days). Additional discrepancies could stem from dif- 
ferent characterizations of completeness, reliance on pho- 
tometric versus spectroscopic measurements of R^^ and 
magnitude-limited, rather than noise-limited, samples. 

9. DISCUSSION 

9.1. TERRA 



We implement in this work a new pipeline for the de- 
tection of transiting planets in Kepler photometry and 
apply it to a sample of 12,000 G and K-type dwarfs stars 
chosen to be among the most photometrically quiet of the 
Kepler target stars. These low noise stars offer the best 
chance for the detection of small. Earth-size planets in 
the Kepler field and will one day be among the stars from 
which ?70 — the fraction of Sun-like stars bearing Earth- 
size planets in habitable zone orbits — is estimated. In 
this work, we focus on the close-in planets having orbital 
periods of 5-50 days and semi-major axes ^ 0.25 AU. 
Earth-size planets with these characteristics are statisti- 
cally at the margins of detectability with the current ~ 
3 years of photometry in Kepler quarters Q1-Q12. 

Our TERRA pipeline has two key features that enable 
confident measurement of the occurrence of close-in plan- 
ets approaching Earth size. First, TERRA calibrates the 
Kepler photometry and searches for transit signals in- 
dependent of the results from the Kepler Mission's offi- 
cial pipeline. In some cases, TERRA calibration achieves 
superior noise suppression compared to the Pre-search 
Data Conditioning (PDC) m odule of the official Kepler 
pipeline (Petigura fc Marcv 2012). The transit search al- 
gorithm in TERRA is efficient at detecting low SNR tran- 
sits in the calibrated light curves. This algorithm suc- 
cessfully redis c overs 82 of 86 stars bearing planets in 
IBatalha et al.l (|2012l ). Recall that the current version 
of TERRA only detects the highest SNR transit signal in 
each system. Thus, additional planets orbiting known 
hosts are not reported here. We report the occurrence of 
stars having one or more planets, not the mean number 
of planets per star as in H12 and el sewhere. Our p i peline 
also detects 37 planets not found in lBatalha et al.l (|2012l ) 
(19 of which were not in the catalog of the Kepler team as 
of August 8*^, 2012), albeit with the benefit of 6 quarters 
of additional photometry for TERRA to search. 
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Figure 15. As in Figure [Tol red points show 119 TERRA-detected planets. Blue points represent additional planets from Batalha et al.l 
{2012). Most (28 out of 32) of these new candidates are planets in multi-candidate systems where TERRA successfully identifies the higher 
SNR candidate. We apply no completeness correction to these new planets, and we believe this is appropriate for Rp > 1.4 Re- We quote 
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The second crucial feature of TERRA is that we have 
characterized its detection completeness via the injec- 
tion and recovery of synthetic transits in real Kepler 
light curves from the Best 12k sample. This complete- 
ness study is crucial to our occurrence calculations be- 
cause it allows us to statistically correct for incomplete- 
ness variations across the P-Rp plane. While the Kepler 
Project has initiated a completeness study of the offi- 
cial pipeline (Christiansen et alj |2012f ), TERRA is the only 
pipeline for Kepler photometry whose detection com- 
pleteness has been calibrated by injection and recovery 
tests. Prior to TERRA, occurrence calculations required 
one to assume that the Kepler planet detections were 
complete down to some SNR limit, or to estimate com- 
pleteness based on SNR alone without empirical tests of 
the performance of the algorithms in the pipeline. For 
example, H12 made cuts in stellar brightness {Kp < 15) 
and transit SNR (> 10 in a single quarter of photom- 
etry) and restricted their search to planets larger than 
2 Re with orbital periods shorter than 50 days. These 



conservative cuts on the planet and star catalogs were 
driven by the unknown completeness of the official Ke- 
pler pipeline at low SNR. HI 2 applied two statistical 
corrections to convert their distribution of detected plan- 
ets into an occurrence distribution. They corrected for 
non-transiting planets with a geometric a/i?^ correction. 
They also computed the number of stars amenable to 
the detection (at SNR > 10 in a single quarter) of each 
planet and considered only that number of stars in the 
occurrence calculation. H12 had no empirical way to de- 
termine the actual detection efficiency of the algorithms 
in the pipeline. Here, we apply the geometric a/i?^ cor- 
rection and correct for pipeline completeness across the 
P-Rp plane by explicit tests of the TERRA pipeline effi- 
ciency, which naturally incorporates an SNR threshold 
correction as in HI 2. 

9.2. Planet Occurrence 

H12 found that for close-in planets, the planet radius 
function rises steeply from Jupiter size to 2 Re- For 
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Figure 16. Same as Figure [TTI with inclusion of planets in multi 
planet systems. The blue regions represent the additional contri- 
bution to planet occurrence from the Bataiha et al. (2012) planets. 
The addition of these new planets does not change the overall shape 
of the distribution. The dashed line is the power law fit to the 
planet size distribution in H12. The fit agrees qualitatively for Rp 
> 2 Re, but not within errors. We expect the H12 fit to be ~ 25% 
higher than our occurrence measurements since H12 included plan- 
ets with P < 50 days (not P = 5-50 days). Additional discrepancies 
between occurrence in H12 and this work could stem from different 
characterizations of completeness, reliance on photometric versus 
spectroscopic measurements of i?^, and magnitude-limited, rather 
than noise-limited, samples. 

smaller planets of ~ 1-2 Re, occurrence was approxi- 
mately constant in logarithmic Rp bins, but HI 2 were 
skeptical of the result below 2 Re because of unknown 
pipeline completeness and the s mall number statistics 
near 1 Re in the Bor ucki et alJ (|2Qllb ) planet catalog. 
In this work, we strongly confirm the power law rise in 
occurrence from 4 to 2 Re using a superior assessment 
of completeness and nine times more photometry than 
in H12. Using TERRA, we can empirically and confidently 
compute occurrence down to 1 Re- Our key result is 
the plateau of planet occurrence for the size range 1-2.8 
Re for planets having orbital periods 5-50 days around 
Sun-like stars. In that size range of 1-2.8 Re, 23% of 
stars have a planet orbiting with periods between 5 and 
50 days. Including the multiple planets within each sys- 
tem, we find 0.28 planets per star within the size range 
1-2.8 Re and with periods between 5-50 days. These 
results apply, of course, to the Kepler field, with its still 
unknown distribution of masses, ages, and metallicities 
in the Galactic disk. 

As shown in Figure [TOl TERRA detects many sub-Earth 
size planets (< 1.0 Re)- These sub-Earths appear in re- 
gions of low completeness, and, provocatively, may repre- 
sent just the tip of the iceberg. A rich population of sub- 
Earths may await discovery given more photometry and 
continued pipeline improvements. With 8 years of total 
photometry in an extended Kepler mission (compared to 
3 years here), the computational machinery of TERRA — 
including its light curve calibration, transit search, and 
completeness calibration — will enable a measurement of 
7^0 for habitable zone orbits. 

9.3. Interpretation 

We are not the first to note the huge population 
of close-in planets smaller o r less massive than Nep- 
tune. Usin^ Doppler surveys, iHoward et al.l (|201Ql ) and 
iMavor et al.l (gOll) showed that the planet mass func- 



tion rises steeply with decreasing mass, at least for close- 
in planets. In Kepler data, the excess of close-in, small 
planets was obvious in th e initial planet catalogs released 
by the Kepler Project (|Borucki et al.1 2011a.b). H12 
characterized the occurrence distribution of these small 
planets as a function of their size, orbital period, and 
host star temperature. These occurrence measurements, 
based on official Kepl er pla n et catalogs, w e re refined and 
extended by Youdi nI (|2QTTI) iTraubl (|2Q12l ). iDong fc Zhul 
([2012|) , [Beauge & Nesvornvl (|201 3V and others. Our con- 
tribution here shows a clear plateau in occurrence in 
the 1-2.8 Re size range and certified by an independent 
search of Kepler photometry using a pipeline calibrated 
by injection and recovery tests. The onset of the plateau 
at ~ 2.8 Re suggests that there is a preferred size scale 
for the for mation of close -in planets. 

HI 2 and YoudinI (|2011f ) noted falling planet occurrence 
for periods shorter than ^ 7 days. We also observe declin- 
ing planet occurrence for short orbital periods, but find 
that the transition occurs closer to ~ 10 days. We con- 
sider our period distribution to be i n qu alitative agree- 
ment with those of HI 2 and lYoudinl (|20 11). Planet for- 
mation and/or migration seems to discourage very close- 
in planets (P ^ 10 days). 

Close-in, small planets are now the most abundant 
planets detected by current transit and Doppler searches, 
yet they are absent from the solar system. The solar sys- 
tem is devoid of planets between 1 and 3.88 Re (Earth 
and Neptune) and planets with periods less than Mer- 
cury's (P < 88.0 days). The formation mechanisms and 
possible subsequent migration of such planet s are hotly 
debat ed. T he population syn t hesis models of llda fc LinI 
(|2010[ ) and iMordasini eTaLl (|2012f ) suggest that they 
form near or beyond the ice line and then migrate qui- 
escently in the protoplanetary disk. These models follow 
the growth and migration of planets over a wide range 
of parameters (from Jupiter mass down to Earth mass 
orbiting at distances out to '^ 10 AU) and they predict 
"deserts" of planet oc currence that are n ot detected. 

Mor e recently, Hansen fc Murravl (12012) and 
Chian g^ Laughhnl (|2012f ) have argued for the in 
situ formation of close-in planets of Neptune size and 
smaller. In these models, close-in rocky planets of a 
few Earth masses form from protoplanetary disks more 
massive than the minimum mass solar nebula. Multiple 
planets per disk form commonly in these models and 
accretion is fast (~ 10^ years) and efficient due to 
the short dynamical timescales of close-in orbits. The 
rocky cores form before the protoplanetary disk has 
dissipated, accreting nebular gas that adds t ypically 
r^ 3% to the mass of the planet (Chiang fc LaughlinI 
l2012f ). But the small amounts of gas can significantly 
swell the radii of thes e otherwise rocky planets. For 
example, I Adams et al.l (;2008.) found that adding a H/He 
gas envelope equivalent to 0.2-20% of the mass of a 
solid 5 Me planet increases the radius 8-110% above 
the gas-free value. 

We find the in situ model plausible because it natu- 
rally explains the large number of close, sub-Neptune- 
size planets, the high rate of pla net multiplicity and 
nearly co-planar and c ircular orbits ( jLissauer et al.|[2QTll : 
iFang fc Margotll2012l ), and does not require tuning of 
planet migration models. Our result of a plateau in the 
planet size distribution for 1-2.8 Re with a sharp falloff 
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in occurrence for larger planets along with decreasing 
occurrence for P ^ 10 days are two significant observed 
properties of planets around Sun-like stars that must be 
reproduced by models that form planets in situ or oth- 
erwise and by associated population synthesis models. 

The in situ model seems supported by the sheer large 
occurrence of sub-Neptune-size planets within 0.25 AU. 
It seems unlikely that all such planets form beyond the 
snow line at ~ 2 AU, which would require inward migra- 
tion to within 0.25 AU, but not all the way into the star. 
Such models of formation beyond the snow line seem 
to require fine tuning of migration and parking mech- 
anisms, as well as the tuning of available water or gas 
beyond 2 AU, while avoiding runaway gas accretion to- 
ward Jupiter masses. Still, in situ formation seems to 
require higher densities than t hose normally assumed in 
a minimum mass solar nebula (iChiang fc Laughlinll2012[ ) 
in order to form the sub-Neptune planets before removal 
of the gas. If this in situva.ode\ is correct, we expect these 
sub-Neptune-size planets to be com posed of rock plus H 
and He, rather than rock plus water (jChiang fc LaughlinI 
[2Q12l) . Thus, a test of the in situ mode of formation in- 
volves spectroscopic measurements of the chemical com- 
position of the close-in sub-Neptunes. 
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Planet candidates identified with TERRA 
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69.883 


1.07 


0.11 


1.35 


0.08 


< 


0.65 




5302 


3.83 


2.45 


2.86 


0.29 


s 


N 


Y 


4049901 


16.291 


65.115 


0.64 


0.23 


1.72 


0.16 


< 


0.68 




5250 


4.48 


0.85 


0.60 


0.21 


PI 


N 


Y 


4548011 


6.284 


66.198 


0.60 


0.06 


2.23 


0.64 


< 


0.49 




5991 


4.30 


1.26 


0.83 


0.09 


s 




N 


4644604 


14.486 


64.550 


2.03 


0.21 


1.49 


0.15 


< 


0.59 




5739 


4.34 


1.04 


2.32 


0.23 


s 


N 


Y 


4770174 


6.096 


67.600 


0.57 


0.20 


3.08 


0.44 


< 


0.75 




6013 


4.44 


1.01 


0.63 


0.22 


PI 


N 


N 


4827723 


7.239 


68.024 


1.63 


0.17 


2.20 


0.37 




0.70 


0.12 


5392 


4.52 


0.87 


1.54 


0.16 


s 


N 


Y 


4914423 


15.965 


75.182 


2.17 


0.22 


3.14 


0.28 




0.71 


0.06 


5904 


4.27 


1.29 


3.05 


0.31 


s 


N 


Y 


4914566 


22.241 


77.736 


0.83 


0.29 


3.76 


0.36 


< 


0.73 




5974 


4.22 


1.31 


1.18 


0.41 


PI 


Y 


N 


5009743 


41.699 


102.563 


1.94 


0.20 


2.87 


0.30 


< 


0.70 




5937 


4.35 


1.09 


2.32 


0.23 


s 


N 


Y 


5042210 


12.147 


64.535 


0.81 


0.08 


3.22 


0.21 


< 


0.56 




6007 


4.27 


1.31 


1.16 


0.12 


s 


N 


Y 


5094751 


6.482 


68.943 


1.69 


0.17 


2.51 


0.31 




0.69 


0.13 


5929 


4.37 


1.10 


2.02 


0.20 


s 


N 


Y 


5096590 


29.610 


70.332 


1.00 


0.35 


2.60 


0.20 


< 


0.67 




5623 


4.63 


0.73 


0.79 


0.28 


PI 




N 


5121511 


30.996 


93.790 


2.35 


0.24 


1.40 


0.06 


< 


0.47 




5217 


4.53 


0.84 


2.15 


0.22 


s 


N 


Y 


5308537 


14.265 


76.217 


0.76 


0.08 


1.92 


0.53 


< 


0.87 




5831 


4.29 


1.26 


1.05 


0.11 


s 




N 


5561278 


20.310 


79.826 


1.16 


0.12 


2.73 


0.15 


< 


0.47 




6161 


4.38 


1.14 


1.45 


0.15 


s 


N 


Y 


5613330 


23.449 


69.031 


1.68 


0.17 


4.44 


0.28 




0.25 


0.22 


6080 


4.20 


1.48 


2.70 


0.27 


s 


N 


Y 


5652893 


14.010 


65.959 


1.15 


0.13 


1.98 


0.49 


< 


0.85 




5150 


4.55 


0.78 


0.97 


0.11 


s 


N 


N 


5702939 


18.398 


77.934 


1.09 


0.38 


2.45 


0.27 


< 


0.54 




5634 


4.47 


0.87 


1.04 


0.36 


PI 


N 


N 


5735762 


9.674 


68.009 


2.73 


0.27 


1.61 


0.06 


< 


0.47 




5195 


4.53 


0.84 


2.51 


0.25 


s 


N 


Y 


5866724 


5.860 


65.040 


1.65 


0.17 


2.01 


0.05 




0.22 


0.14 


6109 


4.29 


1.31 


2.37 


0.24 


s 


N 


Y 


5959719 


6.738 


66.327 


0.99 


0.10 


1.30 


0.21 


< 


0.76 




5166 


4.56 


0.77 


0.83 


0.09 


s 


N 


Y 


6071903 


24.308 


87.054 


2.28 


0.23 


1.29 


0.05 


< 


0.27 




5296 


4.55 


0.83 


2.06 


0.21 


s 


N 


Y 


6197215 


10.613 


68.691 


1.23 


0.13 


0.48 


0.06 


< 


0.80 




5933 


4.39 


1.09 


1.46 


0.15 


s 


N 


N 


6289257 


19.675 


69.903 


1.33 


0.13 


1.96 


0.21 




0.35 


0.27 


6023 


4.36 


1.11 


1.61 


0.16 


s 


N 


Y 


6291837 


35.596 


84.942 


2.45 


0.25 


3.18 


0.42 


< 


0.71 




6165 


4.38 


1.14 


3.05 


0.31 


s 


N 


Y 


6356692 


11.392 


74.555 


0.66 


0.23 


2.74 


0.31 


< 


0.73 




5420 


4.03 


1.64 


1.19 


0.42 


PI 


N 


N 


6365156 


10.214 


73.050 


1.57 


0.16 


2.83 


0.09 


< 


0.42 




5852 


4.28 


1.25 


2.13 


0.21 


s 


N 


Y 


6442340 


13.137 


76.973 


1.38 


0.14 


2.41 


0.13 


< 


0.55 




5764 


4.38 


1.08 


1.63 


0.16 


s 


N 


Y 


6521045 


12.816 


68.772 


1.37 


0.14 


3.42 


0.34 




0.40 


0.19 


5874 


4.29 


1.24 


1.85 


0.19 


s 


N 


Y 


6523351 


6.067 


69.128 


0.72 


0.25 


1.02 


0.20 


< 


0.89 




5489 


4.15 


1.35 


1.06 


0.37 


PI 




N 


6605493 


9.310 


69.379 


0.99 


0.10 


1.80 


0.16 


< 


0.62 




5805 


4.36 


1.12 


1.20 


0.12 


s 


N 


Y 


6607357 


7.700 


67.390 


0.85 


0.30 


2.72 


0.45 




0.55 


0.22 


5592 


4.51 


0.91 


0.84 


0.29 


PI 


N 


N 


6707835 


22.248 


84.881 


2.25 


0.23 


1.97 


0.07 


< 


0.43 




5619 


4.44 


0.99 


2.42 


0.24 


s 


N 


Y 


6716545 


13.910 


75.855 


0.84 


0.30 


3.29 


0.68 


< 


0.85 




6044 


4.30 


1.12 


1.03 


0.36 


PI 


N 


N 


6803202 


21.061 


76.591 


1.58 


0.16 


2.64 


0.08 


< 


0.36 




5719 


4.40 


1.03 


1.77 


0.18 


s 


N 


Y 


6851425 


11.120 


72.744 


2.29 


0.23 


1.71 


0.12 


< 


0.62 




5071 


4.59 


0.72 


1.80 


0.18 


s 


N 


Y 


6922710 


23.127 


78.663 


1.01 


0.11 


2.57 


0.47 


< 


0.78 




5929 


4.40 


1.07 


1.18 


0.12 


s 


N 


Y 


7021534 


9.066 


68.147 


1.42 


0.50 


0.81 


0.05 


< 


0.55 




5848 


4.55 


0.87 


1.35 


0.47 


PI 


Y 


N 


7033671 


9.490 


66.961 


1.48 


0.15 


1.81 


0.09 


< 


0.46 




5679 


4.29 


1.26 


2.03 


0.20 


s 


N 


Y 


7211221 


5.621 


69.903 


1.21 


0.12 


1.23 


0.11 


< 


0.53 




5634 


4.44 


0.93 


1.23 


0.12 


s 


N 


Y 


7219825 


17.233 


68.084 


2.06 


0.21 


2.27 


0.13 


< 


0.22 




6089 


4.36 


1.13 


2.54 


0.25 


s 


N 


Y 


7345248 


5.665 


69.338 


0.63 


0.22 


2.81 


0.43 


< 


0.61 




5656 


4.27 


1.19 


0.81 


0.29 


P2 




N 


7419318 


18.736 


73.153 


2.13 


0.22 


1.47 


0.11 


< 


0.61 




5187 


4.54 


0.81 


1.89 


0.19 


s 


N 


Y 


7466863 


11.971 


68.173 


1.83 


0.64 


0.98 


0.03 


< 


0.39 




6035 


4.12 


1.45 


2.89 


1.01 


PI 


Y 


N 


7582689 


11.921 


70.475 


0.75 


0.26 


5.84 


1.24 


< 


0.95 




6022 


4.04 


1.64 


1.34 


0.47 


PI 




N 


7668663 


6.498 


69.012 


1.44 


0.15 


1.85 


0.34 




0.79 


0.09 


5725 


4.33 


1.14 


1.79 


0.18 


s 


N 


Y 


7700622 


35.585 


86.426 


2.83 


0.28 


1.99 


0.10 


< 


0.46 




4787 


4.62 


0.69 


2.13 


0.21 


s 


N 


Y 


7762723 


9.887 


72.539 


0.73 


0.26 


2.21 


0.39 


< 


0.65 




5501 


4.58 


0.80 


0.64 


0.22 


PI 


Y 


N 


7810483 


29.921 


79.024 


1.71 


0.60 


1.58 


0.11 


< 


0.58 




5893 


4.52 


0.91 


1.70 


0.59 


PI 


Y 


N 


7906739 


7.015 


69.903 


0.90 


0.31 


2.43 


0.32 


< 


0.68 




5652 


4.53 


0.81 


0.80 


0.28 


PI 


Y 


N 


7906892 


8.849 


72.797 


0.58 


0.07 


8.03 


3.69 




0.89 


0.12 


6095 


4.35 


1.14 


0.72 


0.08 


s 




N 


7918652 


11.456 


69.316 


0.79 


0.28 


2.45 


0.47 


< 


0.81 




5809 


4.25 


1.19 


1.02 


0.36 


PI 


N 


N 


8008067 


15.771 


70.584 


2.22 


0.22 


3.32 


0.31 




0.67 


0.08 


5594 


4.37 


1.10 


2.66 


0.27 


s 


N 


Y 


8009496 


38.476 


83.567 


1.88 


0.66 


1.72 


0.14 


< 


0.70 




5833 


4.54 


0.85 


1.74 


0.61 


PI 


Y 


N 


8073705 


10.601 


65.967 


0.75 


0.08 


2.29 


0.36 


< 


0.80 




6086 


4.36 


1.12 


0.91 


0.10 


s 




N 


8077137 


15.090 


78.772 


0.78 


0.08 


2.47 


0.65 


< 


0.85 




6179 


4.37 


1.16 


0.99 


0.10 


s 


N 


Y 


8081187 


37.323 


85.828 


1.67 


0.59 


3.58 


0.39 


< 


0.62 




6030 


4.55 


0.86 


1.57 


0.55 


PI 


Y 


N 


8087812 


27.211 


65.360 


1.01 


0.10 


5.01 


0.57 


< 


0.67 




5985 


4.17 


1.41 


1.55 


0.16 


s 




N 


8242434 


44.964 


77.565 


2.58 


0.26 


2.58 


0.15 




0.31 


0.22 


4692 


4.63 


0.68 


1.92 


0.19 


s 


N 


Y 


8280511 


10.435 


67.826 


1.37 


0.14 


1.70 


0.12 


< 


0.58 




5522 


4.45 


0.91 


1.36 


0.14 


s 




N 


8323753 


6.714 


67.308 


1.97 


0.69 


1.33 


0.05 


< 


0.52 




5817 


4.23 


1.26 


2.71 


0.95 


PI 


Y 


N 


8349582 


11.523 


64.959 


2.14 


0.22 


2.41 


0.21 




0.61 


0.09 


5668 


4.23 


1.35 


3.15 


0.32 


s 


N 


Y 


8429668 


5.007 


67.696 


0.73 


0.26 


1.66 


0.47 




0.62 


0.29 


5034 


4.63 


0.65 


0.52 


0.18 


P2 




N 
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Light Curve Fit 










Stellar Parameters 
Teff log^ R^ 


Rp 


criRp) 


source^ 


FP 




kic 


P 


to" 




-(If) 


r 


a{T) 




b 


a{b) 


B12 




(d) 


(d) 


(%) 




(hrs) 










(K) 


(cgs) 


(^o) 


(Re) 










8480285 


29.667 


92.687 


2.43 


0.25 


4.72 


0.42 




0.48 


0.14 


5960 


4.36 


1.09 


2.89 


0.29 


S 


N 


Y 


8494617 


22.923 


66.640 


1.06 


0.11 


3.72 


0.36 


< 


0.50 




5905 


4.36 


1.10 


1.27 


0.13 


s 


N 


Y 


8560804 


31.976 


66.803 


0.99 


0.35 


4.96 


0.48 


< 


0.65 




5878 


4.44 


1.01 


1.09 


0.38 


PI 


N 


N 


8611832 


22.597 


73.431 


1.08 


0.11 


3.15 


0.17 


< 


0.58 




5577 


4.37 


0.99 


1.17 


0.12 


s 


N 


Y 


8628758 


14.374 


71.204 


1.99 


0.20 


4.76 


0.68 


< 


0.91 




5773 


4.40 


1.04 


2.26 


0.23 


s 


N 


Y 


8631504 


14.820 


66.409 


1.20 


0.12 


2.02 


0.25 


< 


0.63 




4828 


4.60 


0.72 


0.95 


0.10 


s 


N 


Y 


8644365 


19.917 


72.354 


1.07 


0.11 


3.37 


0.55 


< 


0.85 




6054 


4.39 


1.12 


1.31 


0.13 


s 




N 


8804455 


7.597 


64.958 


1.14 


0.12 


2.55 


0.47 




0.74 


0.11 


5715 


4.38 


1.07 


1.33 


0.14 


s 


N 


Y 


8805348 


29.907 


78.383 


2.36 


0.24 


3.20 


0.39 




0.66 


0.12 


5739 


4.34 


1.05 


2.69 


0.27 


s 


N 


Y 


8822366 


30.864 


65.012 


1.41 


0.14 


4.16 


0.20 


< 


0.62 




6089 


4.35 


1.14 


1.76 


0.18 


s 


N 


Y 


8827575 


10.129 


68.997 


0.84 


0.30 


1.94 


0.18 


< 


0.63 




5284 


4.45 


0.89 


0.82 


0.29 


PI 




N 


8866102 


17.834 


114.225 


1.66 


0.17 


2.28 


0.06 


< 


0.39 




6178 


4.37 


1.15 


2.08 


0.21 


s 


N 


Y 


8962094 


30.865 


75.052 


2.09 


0.21 


1.37 


0.10 


< 


0.61 




5739 


4.34 


1.04 


2.38 


0.24 


s 


N 


Y 


8972058 


8.991 


69.746 


2.01 


0.20 


2.10 


0.11 


< 


0.52 




5979 


4.38 


1.09 


2.39 


0.24 


s 


N 


Y 


9006186 


5.453 


67.152 


0.85 


0.09 


1.07 


0.08 


< 


0.61 




5404 


4.53 


0.87 


0.81 


0.08 


s 


N 


Y 


9086251 


6.892 


64.632 


0.87 


0.09 


0.93 


0.10 


< 


0.68 




6044 


4.22 


1.45 


1.38 


0.14 


s 


N 


Y 


9139084 


5.836 


67.853 


2.07 


0.21 


1.06 


0.03 


< 


0.41 




5411 


4.53 


0.85 


1.92 


0.19 


s 


N 


Y 


9226339 


21.461 


65.230 


1.10 


0.11 


1.66 


0.16 


< 


0.70 




5807 


4.28 


1.25 


1.50 


0.15 


s 




N 


9288237 


7.491 


68.329 


0.52 


0.18 


3.01 


0.73 


< 


0.84 




5946 


4.44 


0.96 


0.54 


0.19 


PI 




N 


9491832 


49.565 


103.693 


1.14 


0.12 


6.43 


1.49 




0.74 


0.25 


5821 


4.15 


1.57 


1.95 


0.21 


s 




N 


9549648 


5.992 


69.883 


1.46 


0.15 


1.80 


0.45 




0.89 


0.06 


6165 


4.38 


1.14 


1.82 


0.19 


s 


N 


Y 


9704384 


5.509 


65.285 


1.35 


0.14 


1.85 


0.26 


< 


0.75 




5448 


4.50 


0.91 


1.35 


0.14 


s 


N 


Y 


9716028 


17.373 


71.258 


0.82 


0.08 


2.46 


0.27 


< 


0.73 




6119 


4.37 


1.15 


1.03 


0.11 


s 




N 


9717943 


6.110 


69.903 


0.72 


0.08 


1.80 


0.87 




0.79 


0.20 


5968 


4.30 


1.27 


1.00 


0.11 


s 


N 


Y 


9886361 


7.031 


67.464 


0.88 


0.09 


3.04 


0.21 


< 


0.42 




6090 


4.39 


1.13 


1.08 


0.11 


s 


N 


N 


10055126 


9.176 


71.722 


1.33 


0.13 


2.43 


0.24 


< 


0.48 




5905 


4.36 


1.10 


1.59 


0.16 


s 


N 


Y 


10130039 


12.758 


66.961 


1.19 


0.12 


2.19 


0.07 


< 


0.38 




5828 


4.42 


1.01 


1.31 


0.13 


s 


N 


Y 


10136549 


9.693 


65.809 


1.14 


0.12 


3.44 


0.44 




0.57 


0.20 


5684 


4.13 


1.59 


1.98 


0.20 


s 


N 


Y 


10212441 


15.044 


66.211 


0.95 


0.10 


2.97 


0.24 


< 


0.57 




5939 


4.35 


1.09 


1.13 


0.11 


s 


N 


Y 


10593535 


20.925 


67.900 


0.93 


0.10 


3.63 


0.59 


< 


0.81 




5822 


4.28 


1.25 


1.27 


0.13 


s 




N 


10722485 


7.849 


67.907 


0.87 


0.09 


2.61 


0.51 


< 


0.87 




5682 


4.36 


1.03 


0.98 


0.10 


s 




N 


10917433 


6.912 


65.190 


0.51 


0.05 


2.03 


0.57 


< 


0.87 




5680 


4.33 


1.12 


0.62 


0.06 


s 




N 


11086270 


31.720 


75.821 


1.90 


0.19 


3.35 


0.40 




0.70 


0.14 


5960 


4.37 


1.08 


2.24 


0.22 


s 


N 


Y 


11121752 


7.630 


70.087 


1.00 


0.10 


1.56 


0.23 


< 


0.75 




6045 


4.36 


1.12 


1.22 


0.12 


s 


N 


Y 


11133306 


41.746 


101.657 


1.90 


0.19 


2.32 


0.15 




0.27 


0.22 


5953 


4.37 


1.10 


2.27 


0.23 


s 


N 


Y 


11241912 


14.427 


71.857 


0.95 


0.10 


2.40 


0.18 


< 


0.59 




5931 


4.40 


1.07 


1.10 


0.11 


s 




N 


11250587 


7.257 


67.028 


1.95 


0.20 


2.41 


0.05 


< 


0.40 




5853 


4.18 


1.49 


3.18 


0.32 


s 


N 


Y 


11253711 


17.791 


82.259 


1.86 


0.65 


1.29 


0.09 


< 


0.57 




5816 


4.48 


0.95 


1.92 


0.67 


PI 


N 


Y 


11295426 


5.399 


69.065 


1.89 


0.19 


2.78 


0.20 




0.80 


0.05 


5793 


4.25 


1.30 


2.68 


0.27 


s 


N 


Y 


11402995 


10.061 


71.959 


2.00 


0.20 


2.42 


0.19 


< 


0.62 




5709 


4.30 


1.18 


2.56 


0.26 


s 


N 


Y 


11554435 


9.434 


73.119 


5.69 


0.57 


1.29 


0.02 


< 


0.43 




5536 


4.52 


0.90 


5.60 


0.56 


s 


N 


Y 


11560897 


35.968 


71.667 


1.46 


0.15 


1.47 


0.11 


< 


0.45 




5832 


4.27 


1.28 


2.03 


0.21 


s 


N 


Y 


11612280 


9.406 


70.699 


0.80 


0.08 


3.02 


0.54 


< 


0.77 




5857 


4.24 


1.32 


1.15 


0.12 


s 




N 


11771430 


40.031 


81.863 


1.39 


0.14 


2.25 


0.14 


< 


0.57 




5850 


4.23 


1.38 


2.10 


0.21 


s 


N 


Y 


11774991 


37.815 


74.112 


1.45 


0.15 


2.37 


0.21 


< 


0.65 




4710 


4.62 


0.70 


1.10 


0.11 


s 


N 


Y 


12254909 


5.350 


66.872 


0.84 


0.29 


2.26 


0.14 


< 


0.60 




5987 


4.46 


0.96 


0.88 


0.31 


PI 


N 


Y 


12301181 


6.147 


67.867 


1.04 


0.11 


1.45 


0.07 


< 


0.40 




4997 


4.60 


0.74 


0.84 


0.09 


s 


N 


Y 


12416661 


8.053 


67.968 


0.63 


0.22 


2.93 


0.37 


< 


0.59 




6091 


4.12 


1.47 


1.02 


0.36 


PI 




N 


12454461 


7.467 


69.392 


0.84 


0.09 


1.83 


0.20 


< 


0.71 




6048 


4.31 


1.27 


1.16 


0.12 


s 


N 


Y 


12737015 


24.669 


69.616 


1.05 


0.11 


4.94 


0.49 


< 


0.69 




6045 


4.15 


1.60 


1.84 


0.19 


s 




N 


Note. 


Orbital 


period, P 


; time of transit center, ^o; planet-to-star radius ratic 


>, Rp/R 


^; the time for the planet to travel 


R^ during 


transit. 


r; and tra 


nsit impact parameter, b are all determined from the[ 
; from SpecMatch. If SpecMatch parameters do not exist 


Mandel & Aeol" r200 


i2^ light curve fit. By default, stel 
:en from the corrected KIC values, 


lar parameters R^, Tgff, 
described in Section|4] ' 


and 


log g come 


,, paramet 


ers are tak 


The 


FP column lists whether a candidate 


was designated a 


false positive by 1 


the Kepler 


^ team (' 


'Y'-yes, 


'N' no. 


' 'no designation) 


. The B12 


column 


lists whether a candidate was 


present 




.etaI.M2012|). 

























Time of transit center (BJD-2454900). 

Source of stellar parameters: 'S' SpecMatch-derived parameters using Keck HIRES spectra, 'PI' photometrically-derived parameters from 
iBatalha et al.l (120121) , 'P2' photometrically-derived parameters computed by the authors. See Section^ for more details. 

Table 3 

New candidates identified with TERRA 



Light Curve Fit 



Stellar Parameters 



kic 



P 

(d) 



to" 
(d) 



if -(if) 
(%) 



(hrs) 



a{r) 



a{b) 



Teff 

(K) 



logs' 

(cgs) 



R^ 

(Rq) 



Rp 

(Re) 



a{Rp) source^ 



3120904 42.915 72. i 



1.15 



0.12 2.96 0.36 < 0.66 



6151 4.31 1.26 1.59 



0.16 
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Table 3 — Continued 



Light Curve Fit 



Stellar Parameters 



kic 


P 


to" 




-(If) 


r 


a{r) 




b a{b) Teff 


logs' 


R^ 


Rp 


criRp) 


source^ 




(d) 


(d) 


(%) 




(hrs) 








(K) 


(cgs) 


(Rq) 


(Re) 






3545135 


8.483 


65.973 


0.81 


0.08 


1.44 


0.10 


<^ 


0.57 


5794 


4.40 


1.02 


0.90 


0.09 


S 


3852655 


11.629 


65.817 


0.85 


0.30 


1.49 


0.23 


< 


0.79 


5822 


4.36 


1.05 


0.97 


0.34 


PI 


4548011 


6.284 


66.198 


0.60 


0.06 


2.23 


0.64 


< 


0.49 


5991 


4.30 


1.26 


0.83 


0.09 


s 


4770174 


6.096 


67.600 


0.57 


0.20 


3.08 


0.44 


< 


0.75 


6013 


4.44 


1.01 


0.63 


0.22 


PI 


5096590 


29.610 


70.332 


1.00 


0.35 


2.60 


0.20 


< 


0.67 


5623 


4.63 


0.73 


0.79 


0.28 


PI 


5308537 


14.265 


76.217 


0.76 


0.08 


1.92 


0.53 


< 


0.87 


5831 


4.29 


1.26 


1.05 


0.11 


s 


5652893 


14.010 


65.959 


1.15 


0.13 


1.98 


0.49 


< 


0.85 


5150 


4.55 


0.78 


0.97 


0.11 


s 


5702939 


18.398 


77.934 


1.09 


0.38 


2.45 


0.27 


< 


0.54 


5634 


4.47 


0.87 


1.04 


0.36 


PI 


6197215 


10.613 


68.691 


1.23 


0.13 


0.48 


0.06 


^ 


0.80 


5933 


4.39 


1.09 


1.46 


0.15 


s 


6356692 


11.392 


74.555 


0.66 


0.23 


2.74 


0.31 


< 


0.73 


5420 


4.03 


1.64 


1.19 


0.42 


PI 


6523351 


6.067 


69.128 


0.72 


0.25 


1.02 


0.20 


< 


0.89 


5489 


4.15 


1.35 


1.06 


0.37 


PI 


6607357 


7.700 


67.390 


0.85 


0.30 


2.72 


0.45 




0.55 0.22 5592 


4.51 


0.91 


0.84 


0.29 


PI 


6716545 


13.910 


75.855 


0.84 


0.30 


3.29 


0.68 


< 


0.85 


6044 


4.30 


1.12 


1.03 


0.36 


PI 


7345248 


5.665 


69.338 


0.63 


0.22 


2.81 


0.43 


< 


0.61 


5656 


4.27 


1.19 


0.81 


0.29 


P2 


7582689 


11.921 


70.475 


0.75 


0.26 


5.84 


1.24 


^ 


0.95 


6022 


4.04 


1.64 


1.34 


0.47 


PI 


7906892 


8.849 


72.797 


0.58 


0.07 


8.03 


3.69 




0.89 0.12 6095 


4.35 


1.14 


0.72 


0.08 


s 


7918652 


11.456 


69.316 


0.79 


0.28 


2.45 


0.47 


< 


0.81 


5809 


4.25 


1.19 


1.02 


0.36 


PI 


8073705 


10.601 


65.967 


0.75 


0.08 


2.29 


0.36 


^ 


0.80 


6086 


4.36 


1.12 


0.91 


0.10 


s 


8087812 


27.211 


65.360 


1.01 


0.10 


5.01 


0.57 


< 


0.67 


5985 


4.17 


1.41 


1.55 


0.16 


s 


8280511 


10.435 


67.826 


1.37 


0.14 


1.70 


0.12 


< 


0.58 


5522 


4.45 


0.91 


1.36 


0.14 


s 


8429668 


5.007 


67.696 


0.73 


0.26 


1.66 


0.47 




0.62 0.29 5034 


4.63 


0.65 


0.52 


0.18 


P2 


8560804 


31.976 


66.803 


0.99 


0.35 


4.96 


0.48 


< 


0.65 


5878 


4.44 


1.01 


1.09 


0.38 


PI 


8644365 


19.917 


72.354 


1.07 


0.11 


3.37 


0.55 


< 


0.85 


6054 


4.39 


1.12 


1.31 


0.13 


s 


8827575 


10.129 


68.997 


0.84 


0.30 


1.94 


0.18 


< 


0.63 


5284 


4.45 


0.89 


0.82 


0.29 


PI 


9226339 


21.461 


65.230 


1.10 


0.11 


1.66 


0.16 


< 


0.70 


5807 


4.28 


1.25 


1.50 


0.15 


s 


9288237 


7.491 


68.329 


0.52 


0.18 


3.01 


0.73 


^ 


0.84 


5946 


4.44 


0.96 


0.54 


0.19 


PI 


9491832 


49.565 


103.693 


1.14 


0.12 


6.43 


1.49 




0.74 0.25 5821 


4.15 


1.57 


1.95 


0.21 


s 


9716028 


17.373 


71.258 


0.82 


0.08 


2.46 


0.27 


< 


0.73 


6119 


4.37 


1.15 


1.03 


0.11 


s 


9886361 


7.031 


67.464 


0.88 


0.09 


3.04 


0.21 


^ 


0.42 


6090 


4.39 


1.13 


1.08 


0.11 


s 


10593535 


20.925 


67.900 


0.93 


0.10 


3.63 


0.59 


< 


0.81 


5822 


4.28 


1.25 


1.27 


0.13 


s 


10722485 


7.849 


67.907 


0.87 


0.09 


2.61 


0.51 


< 


0.87 


5682 


4.36 


1.03 


0.98 


0.10 


s 


10917433 


6.912 


65.190 


0.51 


0.05 


2.03 


0.57 


^ 


0.87 


5680 


4.33 


1.12 


0.62 


0.06 


s 


11241912 


14.427 


71.857 


0.95 


0.10 


2.40 


0.18 


< 


0.59 


5931 


4.40 


1.07 


1.10 


0.11 


s 


11612280 


9.406 


70.699 


0.80 


0.08 


3.02 


0.54 


< 


0.77 


5857 


4.24 


1.32 


1.15 


0.12 


s 


12416661 


8.053 


67.968 


0.63 


0.22 


2.93 


0.37 


< 


0.59 


6091 


4.12 


1.47 


1.02 


0.36 


PI 


12737015 


24.669 


69.616 


1.05 


0.11 


4.94 


0.49 


< 


0.69 


6045 


4.15 


1.60 


1.84 


0.19 


s 



Note. — The 37 TERRA candidates not in lBatalha et all (120121) and not listed as false positives by the Kepler team. The column definitions are 
the same as in Table [2] 



Table 4 

Union of IBatalha et al.l (|2012l ') and TERRA planet candidate catalogs 



Batalha 



TERRA 



kic 


KOI 


P 


Rp 


P 


Rp 


2142522 








13.32 


1.11 


2307415 


2053.01 


13.12 


1.65 


13.12 


1.57 


2441495 


166.01 


12.49 


2.70 


12.49 


1.99 


2444412 


103.01 


14.91 


2.97 


14.91 


3.41 


2571238 


84.01 


9.29 


2.53 


9.29 


2.82 


2853446 


1118.01 


7.37 


2.52 


7.37 


1.65 


3098810 


1878.01 


40.81 


3.24 


40.81 


2.75 


3120904 








42.91 


1.59 


3342794 


2278.01 


14.17 


1.99 


14.17 


1.68 


3442055 


1218.01 


29.62 


2.22 


29.62 


1.72 


3531558 


118.01 


24.99 


1.41 


24.99 


1.85 


3545135 








8.48 


0.90 


3835670 


149.01 


14.56 


5.50 


14.56 


4.93 


3839488 


1216.01 


11.13 


1.57 


11.13 


1.64 


3852655 








11.63 


0.97 


3942670 


392.02 


12.61 


1.33 






3942670 


392.01 


33.42 


2.27 


33.42 


2.13 


4043190 


1220.01 


6.40 


1.95 


6.40 


2.86 


4049901 


2295.01 


16.29 


0.63 


16.29 


0.60 


4548011 








6.28 


0.83 


4644604 


628.01 


14.49 


1.87 


14.49 


2.32 


4770174 








6.10 


0.63 


4827723 


632.01 


7.24 


1.46 


7.24 


1.54 
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Table 4 — 


Continued 






Batalha 




TERRA 


kic 


KOI 


P 


Rp 


P 


Rp 


4914423 


108.01 


15.97 


2.94 


15.97 


3.05 


4914566 








22.24 


1.18 


5009743 


1609.01 


41.70 


2.34 


41.70 


2.32 


5042210 


2462.01 


12.15 


1.37 


12.15 


1.16 


5094751 


123.01 


6.48 


2.64 


6.48 


2.02 


5094751 


123.02 


21.22 


2.71 






5096590 








29.61 


0.79 


5121511 


640.01 


31.00 


2.43 


31.00 


2.15 


5308537 








14.27 


1.05 


5446285 


142.01 


10.92 


4.08 






5561278 


1621.01 


20.31 


2.48 


20.31 


1.45 


5613330 


649.01 


23.45 


2.31 


23.45 


2.70 


5652893 








14.01 


0.97 


5702939 








18.40 


1.04 


5735762 


148.02 


9.67 


3.14 


9.67 


2.51 


5735762 


148.03 


42.90 


2.35 






5866724 


85.01 


5.86 


2.35 


5.86 


2.37 


5866724 


85.03 


8.13 


1.41 






5959719 


2498.01 


6.74 


0.78 


6.74 


0.83 


6061773 


2001.01 


8.28 


2.39 






6071903 


306.01 


24.31 


2.28 


24.31 


2.06 


6197215 








10.61 


1.46 


6289257 


307.02 


5.21 


1.15 






6289257 


307.01 


19.67 


1.80 


19.67 


1.61 


6291837 


308.01 


35.60 


3.15 


35.60 


3.05 


6356692 








11.39 


1.19 


6365156 


662.01 


10.21 


2.05 


10.21 


2.13 


6442340 


664.02 


7.78 


1.19 






6442340 


664.01 


13.14 


1.83 


13.14 


1.63 


6442340 


664.03 


23.44 


1.19 






6521045 


41.02 


6.89 


1.23 






6521045 


41.01 


12.82 


2.08 


12.82 


1.85 


6521045 


41.03 


35.33 


1.40 






6523351 








6.07 


1.06 


6605493 


2559.01 


9.31 


0.99 


9.31 


1.20 


6607357 








7.70 


0.84 


6678383 


111.01 


11.43 


2.14 






6678383 


111.02 


23.67 


2.05 






6707835 


666.01 


22.25 


2.56 


22.25 


2.42 


6716545 








13.91 


1.03 


6803202 


177.01 


21.06 


1.84 


21.06 


1.77 


6850504 


70.04 


6.10 


0.91 






6850504 


70.01 


10.85 


3.09 






6850504 


70.05 


19.58 


1.02 






6851425 


163.01 


11.12 


2.27 


11.12 


1.80 


6922710 


2087.01 


23.13 


1.54 


23.13 


1.18 


7021534 








9.07 


1.35 


7033671 


670.01 


9.49 


1.92 


9.49 


2.03 


7211221 


1379.01 


5.62 


1.06 


5.62 


1.23 


7219825 


238.01 


17.23 


2.40 


17.23 


2.54 


7219825 


238.02 


26.69 


1.35 






7345248 








5.66 


0.81 


7419318 


313.02 


8.44 


1.61 






7419318 


313.01 


18.74 


2.20 


18.74 


1.89 


7466863 








11.97 


2.89 


7582689 








11.92 


1.34 


7668663 


1898.01 


6.50 


1.50 


6.50 


1.79 


7700622 


315.01 


35.59 


2.14 


35.59 


2.13 


7762723 








9.89 


0.64 


7810483 








29.92 


1.70 


7906739 








7.01 


0.80 


7906892 








8.85 


0.72 


7918652 








11.46 


1.02 


8008067 


316.01 


15.77 


2.72 


15.77 


2.66 


8009496 








38.48 


1.74 


8073705 








10.60 


0.91 


8077137 


274.01 


15.09 


1.12 


15.09 


0.99 


8077137 


274.02 


22.80 


1.13 






8081187 








37.32 


1.57 


8087812 








27.21 


1.55 


8242434 


1726.01 


44.96 


5.25 


44.96 


1.92 


8280511 


1151.01 


5.22 


0.84 






8280511 


1151.02 


7.41 


0.97 






8280511 








10.44 


1.36 


8323753 








6.71 


2.71 
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Table 4 — Continued 





Batalha 




TERRA 


kic 


KOI 


P 


Rp 


P 


Rp 


8349582 


122.01 


11.52 


2.78 


11.52 


3.15 


8429668 








5.01 


0.52 


8480285 


691.02 


16.23 


1.25 






8480285 


691.01 


29.67 


2.92 


29.67 


2.89 


8494617 


2389.01 


22.92 


1.45 


22.92 


1.27 


8554498 


5.02 


7.05 


0.66 






8560804 








31.98 


1.09 


8611832 


2414.01 


22.60 


1.03 


22.60 


1.17 


8611832 


2414.02 


45.35 


1.17 






8628758 


1279.02 


9.65 


0.74 






8628758 


1279.01 


14.37 


1.31 


14.37 


2.26 


8631504 


2503.01 


14.82 


2.41 


14.82 


0.95 


8644365 








19.92 


1.31 


8804455 


2159.01 


7.60 


1.01 


7.60 


1.33 


8805348 


695.01 


29.91 


2.51 


29.91 


2.69 


8822366 


1282.01 


30.86 


3.00 


30.86 


1.76 


8827575 








10.13 


0.82 


8866102 


42.01 


17.83 


2.71 


17.83 


2.08 


8962094 


700.02 


9.36 


1.29 






8962094 


700.03 


14.67 


1.29 






8962094 


700.01 


30.86 


2.28 


30.87 


2.38 


8972058 


159.01 


8.99 


2.70 


8.99 


2.39 


9006186 


2169.01 


5.45 


1.02 


5.45 


0.81 


9086251 


2367.01 


6.89 


1.17 


6.89 


1.38 


9139084 


323.01 


5.84 


2.17 


5.84 


1.92 


9226339 








21.46 


1.50 


9288237 








7.49 


0.54 


9471974 


119.01 


49.18 


3.76 






9491832 








49.57 


1.95 


9549648 


1886.01 


5.99 


2.45 


5.99 


1.82 


9704384 


1913.01 


5.51 


1.40 


5.51 


1.35 


9716028 








17.37 


1.03 


9717943 


2273.01 


6.11 


1.02 


6.11 


1.00 


9886361 








7.03 


1.08 


10055126 


1608.01 


9.18 


1.81 


9.18 


1.59 


10055126 


1608.02 


19.74 


1.58 






10130039 


1909.02 


5.47 


1.15 






10130039 


1909.01 


12.76 


1.52 


12.76 


1.31 


10130039 


1909.03 


25.10 


1.63 






10136549 


1929.01 


9.69 


2.00 


9.69 


1.98 


10212441 


2342.01 


15.04 


1.22 


15.04 


1.13 


10593535 








20.92 


1.27 


10722485 








7.85 


0.98 


10917433 








6.91 


0.62 


11086270 


124.01 


12.69 


3.00 






11086270 


124.02 


31.72 


3.58 


31.72 


2.24 


11121752 


2333.02 


7.63 


1.63 


7.63 


1.22 


11133306 


276.01 


41.75 


2.49 
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Note. — All iBatalha et al.l ()2012|) candidates with P = 5 50 days belonging to stars in the Bestl2k sample are included. Candidates are 
considered equal if they bel ong to the sa me star and the periods in each catalog agree to better than 0.01 days. Eighty-two candidates appear in 
both catalogs, 33 appear in Batalha et al. (2012) only, and 47 appear in this work only (although 10 were listed as false positives by the Kepler 
team). Differences in Rp between the two catalogs stem from different values of R^. Most TERRA planet candidates have SpecMatch-derived stellar 
parameters which are more accurate thau iBatalha et al.i (i2Q12i) parameters, which were derived from KIC broadband photometry. 
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Figure 17. Phase-folded photometry for 18 of the 37 TERRA planet candidates, not in^Batalha et al.' (2012), ordered according to size. For 
clarity, we show median photometric measurements in 10 min bins. The red lines are the best-fitting Mandel fc Ag ol (^ 2002 ) model. 



22 



E. Petigura & G. Marcy 



vf^' 



'J^^^* 



*% 



KIC-9716028 
P = 17 Ad 
Rp = 1.03 Re 



% * 



^**\ >* ******* 




t ** 

* + + 



KIC-5702939 
P = 18.4c/ 
Rp = 1.04 Re 






+ *l u'lr RQnQK-^T 



KIC-5308537 
P = 14.3c/ 
/?P = 1.05 Re 






¥ 



KIC-6523351 
P = 6.1d 
Rp = 1.06 /?£ 






+ ;***+ 



KIC-9886361 
P = 7.0c/ 
Pp = 1.08 Re 



^^'^* . 



***** *t*%n^ 



+ + + X** 



+ + 



+ *++ 



KIC-8560804 
P = 32.0c/ 
Pp = 1.09 Re 



^^f^3r?^- i^^?f<:^/fm '^^^^•^W 



+*+ + 



KIC-11241912 
P = 14.4c/ 
Pp = 1.10 Re 



KIC-11612280 
P = 9.4c/ 
Pp = 1.15 Pe 



^**^^ 



KIC-6356692 
P = 11.4c/ 
Rp = 1.19 Pe 






■:* -' *. k^is 






tr**-. 



KIC-10593535 
P = 20.9c/ 
Pp = 1.27 Re 



^Jk^i 



*\aIL 



*/t\*:**f\ _li**.**^* >r^j*i)-^*\- ,; 

++ ^+++^ KIC-8644365 * 



P = 19.9c/ 
Pp = 1.31 Pe 






KIC-7582689 
P = 11.9c/ 
Pp = 1.34 Re 



^ffe^ f?^%^ 



KIC-8280511 
P = 10.4c/ 
Pp = 1.36 Re 







KIC-6197215 
P = 10.6c/ 
Pp = 1.46 Re 




% * ** * * * , 

+ + ++ +^ *i* ** * '* ^ * -^ 

* * * +* :-**^ *^''^+ KIC-8087812 






P = 27.2c/ 
Pp = 1.55 Re 

f * 







KIC-3120904 
P = 42.9c/ 
Pp = 1.59 Re 



*+* 



* * *. ^ * *\. * f\ 

** * V:^* , * *^ *\ KIC-12737015 



12737015 
P = 24.7c/ 
Pp = 1.84 Re 



E 
Q. 
Q. 

O 
O 
CM 



^V. 



*t* 




^** 



#1. 



KIC-9491832 
P = 49.6c/ 
Pp = 1.95 Pe 



12 hours + 

Figure 18. Same as Figure fTTI but showing the remaining 19 of the 37 TERRA planet candidates, not in lBatalha et all (|2012l ). 



